Density Perturbations in the Early Universe 



P. G. Miedema* 
Royal Netherlands Military Academy 
Hogeschoollaan 2, NL-^818 CR Breda, The Netherlands 

W. A. van Leeuwcn 
University of Amsterdam, Institute for Theoretical Physics, 
Valckenierstraat 65, NL-1018 XE Amsterdam, The Netherlands 
(Dated: November 7, 2005) 

We propose a way to construct manifestly gauge independent quantities out of the gauge dependent 
quantities occurring in the linearized Einstein equations. Thereupon, we show that these gauge- 
invariant combinations can be identified with measurable perturbations to the particle and energy 
densities. In order to follow the evolution of these quantities in time, we rewrite the full set of 
linearized Einstein equations as a set of four first order, ordinary differential equations for the four 
physical quantities in question, namely 1. the perturbations to the energy density, 2. the particle 
number density, 3. the divergence of the spatial fluid velocity and 4. the local spatial curvature due 
to a density perturbation. 

Next, we recombine the linearized perturbation equations, which still contain gauge functions, in 
such a way that the resulting equations contain only gauge-invariant quantities; as it should be, all 
gauge functions have disappeared. In this way, we have arrived at a formulation of the equations for 
cosmological perturbations that is manifestly gauge- invariant. This set of linearized and manifestly 
gauge-invariant Einstein equations for the gauge-invariant combinations constitute the main result 
of this paper. 

The set of equations for the new gauge-invariant quantities, however complicated it may seem, 
is more tractable than those of earlier treatments. This lucky circumstance entails that, unlike the 
situation encountered in former treatments, we do not need to have recourse to any approxima- 
tion. Therefore, we obtain results that, sometimes, differ slightly from those found so far. As an 
illustration, we consider the equations for two cases that can be treated analytically, namely the 
radiation-dominated and matter-dominated eras of a flat Friedmann-Lemaitre-Robertson- Walker 
universe. 

In the radiation-dominated era we find, for small-scale perturbations, acoustic waves with an 
increasing amplitude, while standard treatments predict acoustic waves with a decaying amplitude. 
For large-scale perturbations we find exactly the same growth rates as in the standard literature. 

Finally we mention that, when considering the non-relativistic limit of the linearized Einstein 
equations we find — not surprisingly — the Poisson equation. Unlike what has been done sofar, 
however, we did not need to limit the discussion to a static universe, i.e., we found the Poisson 
equation to be valid not only in a static, but also in an expanding universe. 

PACS numbers: 98.80.Bp, 04.25. Nx, 98.65.Dx, 98.62. Ai, 98.80.Jk 



I. INTRODUCTION 



In order to link the density perturbations at the end of the period of inflation and the observed temperature 
perturbations in the cosmic 2.7 K microwave background radiation, one needs the linearized Einstein equations. 
These equations, which determine the growth of densities in the expanding universe after the era of inflation, arc the 
main subject of this article. Our final results are contained in the Eqs. (6)-(8) with (205) and (207). 

Mathematically, the problem is very simple, in principle. First, all quantities relevant to the problem are split 
up into two parts: 'a background part' and 'a perturbation part'. The background parts are chosen to satisfy the 
Einstein equations for an isotropic universe, i.e., one chooses for the background quantities the Friedmann-Lemaitre- 
Robertson- Walker-solution. Because of the isotropy, the background quantities depend on the time coordinate t only. 
The perturbation parts are supposed to be small compared to their background counterparts, and to depend on the 
space-time coordinate x = (ct, x). The background and perturbations are often referred to as 'zero order' and 'first 
order' quantities, respectively, and we will use this terminology also in this article. After substituting the sum of the 
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zero order and first order parts of all relevant quantities into the Einstein equations, all terms that are products of 
two or more quantities of first order are neglected. This procedure leads, by construction, to a set of linear differential 
equations for the quantities of first order. The solution of this set of linear differential equations is then reduced to a 
standard problem of the theory of ordinary, linear differential equations. 



A. Historical embedding 



The first systematic and extensive study of cosmological perturbations is due to Lifshitz [1, 2, 3]. Almost half a 
century later, Mukhanov, Fcldman and Brandcnbcrgcr, in their 1992 review article [4] entitled 'Theory of Cosmological 
Perturbations' mention or discuss more than 60 articles on the subject, and, thereupon, suggest their own approach 
to the problem. An article of Chung-Pei Ma and Bertschinger [5], entitled 'Cosmological Perturbation Theory in the 
Synchronous and Conformal Newtonian Gauges' critically reviews the existing literature on the subject, including the 
work of Mukhanov et al. 

For a solid overview of the literature we refer to the references by the aforementioned authors [4, 5]. For a nice 
historical overview of gauge theories we refer to O'Raifeartaigh and Straumann [6]. An actual linking of inflationary 
perturbations and observable radiation background anisotropies has been performed by van Tent [7]. 

The fact that so many studies are devoted to a problem that is nothing but obtaining the solution of a set of 
ordinary, linear differential equations is due to the fact that there are several complicating factors, not regarding the 
mathematics involved, but with respect to the physical interpretation of the solutions. In this article we will try and 
explain what has been done so far, and we will define, at the same time, gauge-invariant expressions for the physical 
perturbations to the energy and particle number densities, which differ from gauge-invariant combinations encountered 
in the literature. In the limit of low velocities, but not necessarily in the limit of slow expansion of the universe, our 
first order equations reduce, not surprisingly, to one single equation, the well-known Poisson equation. We thus find 
that the Poisson equation is valid in an expanding universe. This result, which has not been found earlier, is an 
immediate result from our choice of gauge-invariant quantities. See Sec. IX. 

Let us now try and explain in some detail what exactly is the complicating factor in any treatment of the linearized 
Einstein equations. 



B. Origin of the problem 



At the very moment that one has split up a physical quantity into a zero order and a first order part, one introduces 
an ambiguity. To be explicit, let us consider the two quantities that play the leading roles in the theory of cosmological 
perturbations. They are: the energy density of the universe, e(x), and the particle number density of the universe, 
n(x). The linearized Einstein equations contain as known functions the zero order functions s m (t) and n (0) (t), which 
describe the evolution of the background, i.e., they describe the evolution of the unperturbed universe and they obey 
the unperturbed Einstein equations, and as unknown functions the perturbations e {1) (x) and n (1) (x). The latter are 
the solutions to be obtained from the linearized Einstein equations. The subindexes and 1, which indicate the order, 
have been put between round brackets, in order to distinguish them from tensor indices. In all calculations, products 
of a zero order and a first order quantity arc considered to be of first order, and are retained, whereas products of 
first order quantities are neglected. Having said all this, we can say where and how the ambiguity mentioned above 
arises. 

The linearized Einstein equations, which, as stated above already, are the equations which determine the first order 
quantities e {l) {x) and n^ 1} (x), do not fix these quantities uniquely. In fact, it turns out that next to any solution for 
£ ( i) and of the linearized Einstein equations, there exist solutions of the form 

e (1) (x) = e (1) (x) +ip(x)d e m (t), (la) 
n m {x) = n m {x) + ip(x)d n m (t), (lb) 

which also satisfy the linearized Einstein equations. Here the symbol do stands for the derivative with respect to 
x = ct. In the equations (1) the function VK 2 -) is an arbitrary but 'small' function of the space-time coordinate 
x = (ct,x), i.e., we consider ip( x ) to be of first order. We will derive the equations (1) at a later point in this 
article; here it is sufficient to note that, as (la) and (lb) show, the perturbations e (1) and n {1) are fixed by the 
linearized Einstein equations up to terms that are proportional to an arbitrary, small function ip(x), usually called 
a gauge function in this context. Since a physical quantity, i.e., a directly measurable property of a system, may 
not depend on an arbitrary function, the quantities and cannot be interpreted as the real, physical, values 
of the perturbations in the energy density or the particle number density. But, if and n (1) are not the physical 
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perturbations, what are the real perturbations? This is the notorious 'gauge problem' encountered in any treatment 
of cosmological perturbations. Many different answers to this question can be found in the literature, none of which 
is completely satisfactory, a fact which explains the ongoing discussion on this subject. It is our hope that the answer 
given in the present article will turn out to be the definitive one. 

C. Gauge-invariant combinations 

What we will do in this article is leave ip( x ) undetermined, i.e., we do not choose any particular gauge, but we 
eliminate the gauge function ip(x) from the theory altogether. This is not a new idea: many articles have been 
written in which gauge- invariant quantities, i.e., quantities independent of a gauge function are used. What one 
essentially does in any gauge-invariant approach is to try and construct expressions for the perturbations [like we 
do in Eqs. (2)], which are such that ^>(x) disappears from the defining expressions of the physical quantities. These 
quantities, the perturbations in energy density and particle number density in our case, are then shown to obey a 
set of linear equations, not containing the gauge function ip(x) anymore [see Eqs. (6) and (7) with (205) and (207)]. 
These equations follow, by elimination of the gauge dependent quantities in favor of the gauge-invariant ones, in a 
straightforward way from the usual linearized Einstein equations, which did contain il){x). In this way, the theory 
is no longer plagued by the gauge freedom that is inherent to the original equations and their solutions: ip( x ) nas 
disappeared altogether, as it should, not with brute force, but as a natural consequence of the definitions (2) of the 
perturbations to the energy and the particle number densities. 

As noted above, our approach also belongs to the class of gauge-invariant cosmological theories, where terms are 
added to the perturbations in such a way that they become independent of the particular choice of coordinates. 
However, this can be done, in principle, in infinitely many ways, since any linear combination of gauge-invariant 
quantities is gauge-invariant also. Our treatment distinguishes itself from earlier treatments by the fact that our set 
of first order equations reduces to the usual non-relativistic theory in the limit that the three-part of the cosmological 
fluid velocity four-vector t/^ 1 is small compared to the velocity of light. Consequently, our splitting up of the energy 
density and the particle number density in a zero order and a first order part is such that the first order part reduces 
to the non-relativistic expression. In other words, our treatment of perturbations is 'around' the classical, Newtonian 
theory of gravity. By the way, this is not a conditio sine qua non for a treatment of cosmological perturbations to be 
true, valid or even useful for some particular purpose, but it is a very desirable property. 

We will come back to this point in Sec. IX. 

D. Structure of the article 

This article being quite extended, it might be useful to outline shortly its structure. 

In Sec. IE we define, in Eqs. (2), the physical perturbations to the energy density, ef^ and the particle number 

density nfl } that we propose in this article. In Sec. IG, we explain how the gauge problem actually arises in cosmology 
once one chooses to solve the Einstein equations perturbativcly, i.e., by means of a series expansion. Furthermore, we 
show in this section that the quantities ef^ and nf^, defined in Sec. IE, are gauge-invariant. 

In Sec. I H we further motivate our choice of these physical perturbations efl } and nf^ by noting that they reduce to 
the perturbations as they occur in Newton's theory of gravitation, contrarily to what is the case for the perturbations 
defined in earlier treatments. The proof of this statement is postponed to Sec. IX. In Sec. X we apply our perturbation 
theory to the cases of a radiation-dominated and the matter-dominated universe and we show that our approach 
yields in the radiation-dominated era small-scale growing density perturbations, in contrast to the decaying density 
perturbations found in earlier literature. 

In Sec. II we introduce a particular system of coordinates, the well-known synchronous coordinates, and rewrite the 
Einstein equations with respect to this particular coordinates. In Sec. Ill we restrict the problem of obtaining a solution 
of the perturbed Einstein equations to universes for which the zero order solution is homogeneous and isotropic, the 
so-called Fricdmann-Robertson- Walker or Fricdmann-Lcmaitre-Robertson- Walker universes. This section is rather 
technical: we treat zero order and first order quantities in Sees. Ill A and IIIB, zero order Einstein equations and 
conservation laws in Sees. IIIC 1 and IIIC2, and, finally, first order equations and conservation laws in Sees. HID 1 
and HID 2. 

In Sec. IV the first order equations are split up according to their tensorial character, i.e., whether they are a scalar, 
a vector or a second rank tensor. The tensorial, the vectorial and the scalar parts of the perturbations, and the 
equations they should obey, are derived in Sees. IV A, IV B and IV C, respectively. Since, in this article, we are only 
interested in the energy density and the particle number density, which both are scalar quantities, we only need the 
scalar equations. In Sec. V, the linearized Einstein equations for scalar perturbations found in Sec. IV C are rewritten 
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and simplified: in this section we show that the set of conservation laws and constraint equations can be rewritten as 
a set of four first order ordinary differential equations and one algebraic equation. It is found that only three scalars 
play a role in the theory of density perturbations, implying that there is very little choice for the construction of 
gauge-invariant quantities. 

A summary of all relevant equations is given in Sec. VI. Sec. VII is devoted to the perturbations in some ther- 
modynamical quantities related to the energy density and particle number density perturbations. In Sec. VIII the 
perturbation equations derived in Sec. V are rewritten in a manifestly gauge-invariant form. In Sec. IX we consider 
the Newtonian limit and in Sec. X we apply the evolution equations (209) to the two main eras of a flat flrw universe. 
Finally, in Sec. XI we relate our results to results of earlier work. 

In Appendix A some useful thermodynamic relations are collected. Appendix B is devoted to gauge transformations. 
In Appendix C we show that in the non-relativistic limit there remains some gauge freedom, left over from the general 
theory relativity. Details to the derivations in Sec. VIII are given in Appendix D. In Appendix E we derive the 
standard equations using the perturbation equations (202). Finally, Appendix F gives a list of the symbols and 
notations used in the main text. 



E. Gauge-invariant quantities 



In the existing literature on cosmological perturbations, one has attempted to solve the problem that corresponds 
to the non-uniqueness of the perturbations £ (1) and n (1) (1) in two, essentially different, ways. The first way is to 
impose an extra condition on the gauge field ip{ x ) [8, 9, 10, 11, 12]. Another way to get rid of the gauge field ip(x) 
is to take linear combinations of the matter variables and components of the perturbed metric tensor to construct 
gauge- invariant quantities (Bardecn [13], Mukhanov, et al. [4]). The latter approach is generally considered better 
than the one where one fixes a gauge, because it not only leads to quantities that are independent of an undetermined 
function, as should be the case for a physical quantity, but it also does not rely on any particular choice for the gauge 
function, and, therefore, has less arbitrariness. 

In this article we approach the gauge problem in a Bardeen- or Mukhanov-like way. However, instead of adding 
terms containing perturbations of the metric tensor field g^ v , we add a term that is proportional to (minus) the 
divergence of the (normalized) cosmological four velocity U^: 

g i ._ _ d £ (Q) a r2 s 



dpn m ( 
do9 (0 ) 



(2b) 



where 9 (0) and are the background and perturbation part of the covariant four-divergence c U^.^ of the cosmo- 
logical fluid velocity field U^(x). At first sight, these additional terms look quite artificial. However, in Sec. IG we 
will show how they arise very naturally. 

In Sec. IG, we will show that the quantities ef^ and nfl } do not change if we switch from the old coordinates to 
new coordinates x M according to 

x" =x" -^(x), (3) 

where the £ p (x) (ix = 0, . . . , 3) are four arbitrary functions, considered to be of first order, of the old coordinates x M . 
In other words, we will first show that 

4(x) = 4(x), (4a) 
nfl(x) = n*(x), (4b) 

i.e., the perturbations (2) are independent of £ M (x), i.e., gauge- invariant. The combinations (2) are not completely 
obvious indeed, but the proof in Sec. I G that they will be gauge-invariant will take away the mystery completely: 
the above gauge-invariant combinations (2) will then even turn out the 'obvious' gauge-invariant quantities to study. 
Before showing this, we round off the present reasoning. 

If we transform the linearized Einstein equations with e (1) (x) and n ( i)(x) to the new coordinates x, (3), it will turn 
out that we find equations in which only the zero component of the gauge functions ^(x), (/i = 0, 1, 2, 3), occurs. We 
will call it ip(x): 

${x) = £°(x). (5) 
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In cosmology, the gauge function ip(x) is to be treated as a first order quantity, i.e., as a small (or 'infinitesimal') 
change of the coordinates. If we eliminate — with the help of (2) — the quantities £(i)(x) and n {1) (x) from the linearized 
Einstein equations in favor of the gauge-invariant quantities £^{x) and nfl^x), we obtain equations in which the gauge 
function i/)(x) is absent altogether. The disappearance of the gauge function ip{x) is due entirely to the introduction 
of the gauge- invariant quantities ef l} (x) and nf l} (x) rather than £ {1) (x) and n {1) (x). 

The gauge-invariant quantities and nf\ could be the physical perturbations of the energy density and the particle 
number density we are after, since they are independent of the gauge field i/j(x), but this is not true a priori, since 
any linear combination of gauge-invariant quantities is gauge-invariant. Hence, it is not evident that the particular 
combinations proposed in the equations (2) correspond to physical, i.e., measurable quantities. However, as we will 
show in section IX, their non-relativistic limit is just what one would expect. 



F. Final results 



The key to the solution of the gauge problem of cosmology are the equations (173)-(174), which are derived in 
Sec. V. We obtain from these equations (see Sec. VIII) a new set of evolution equations for the gauge-invariant 
perturbations efj, and nfl y The equation for ef^ is a simple second order differential equation of the form: 

efi) + aiefi') + a 2 ef t) = a 3 af^, (6) 

where af^ is given by the function (207). A dot denotes the derivative with respect to x° = ct. The coefficients a\, a-2 
and 03 are complicated combinations of £ (0) , n (0) , the spatial curvature 3 i? (0 ), the Hubble function H and the equation 
of state p m = p(£(o) , n {a) ) and its partial derivatives [see Eq. (205)] . After one has solved this equation for e^ 1 , , one can 
calculate the gauge-invariant perturbation to the particle number density with the help of the formula [see Eq. (189)] 



n S i - (jSi , %)) gi ( 7 \ 

(1) ~ CT(1, + £(0) (1 + U ,) e<1) ' [n 

where w := p (0) /e (0) . The second order differential equation (6) for the first order energy density perturbation efl } 
and the expression (7) for the first order particle density perturbation nf', are the main results of this paper. They 
are obtained without any approximation. 

With their help one can relate the (measurable) fluctuation in the cosmic background temperature [see Eq. (200)] 

T?i = , (8) 

to fluctuations T^ } (to,x) at a time to at the end of the era of inflation. 

In order to get some insight into what these equations imply, we consider two particular cases: the flat flrw universe 
in the radiation-dominated and the matter-dominated eras. As is common usage, we consider the so-called 'contrast 
functions' 5 £ := £f 1 I ) /e (0 ) and 6 n := nf^/n^ [see Eqs. (209)] rather than the fluctuations efl } and nfl } themselves. 
In the radiation-dominated era we find for these contrast functions the system of equations (235) with the exact 
solutions (236) and (244). For large-scale perturbations these solutions are in complete agreement with the solutions 
found in the standard theory. For small-scale perturbations, however, Eqs. (235) yield solutions different from those 
found in the standard theory: we find an oscillating solution with an increasing amplitude, whereas the standard 
result is an oscillating solution with a decaying amplitude. 

In the matter-dominated era we find the equation (265), which is slightly different from the standard equation (279). 
For the density contrast function, S e , we find the solution (270) rather than the solution (280) found in the standard 
literature. 

The purpose of this article is to present a logical and straightforward derivation of the linearized Einstein equations, 
culminating in the equations (6), (7) and (8). It is not our purpose to solve and discuss these equations in detail in 
the various stages of the evolution of the universe, and analyze the possible cosmological consequences that might be 
contained in our equations. 
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G. Gauge-invariant first order perturbations 

We now proceed with the proof that ef^ and nf^ are gauge- invariant. To that end, we start by recalling the defining 
expression for the Lie derivative of an arbitrary tensor field A '" 13 AI ...„ with respect to a vector field £ T (x). It reads 

(•£-£^1) /*•••!/ = A® 1 ^ '^...i/;t£ T 

_ AT— f) fO _ AOt—T £0 

S" 1 )Li"'I/S \T /^'"I^S \T 

+ A a -P T ...„C^ + --- + A a -^... T C;„, (9) 

where the semi-colon denotes the covariant derivative. At the right-hand side, there is a term with a plus sign for 
each lower index and a term with a minus sign for each upper index. Recall also, that the covariant derivative in the 
expression for the Lie derivative may be replaced by an ordinary derivative, since the Lie derivative is, by definition, 
independent of the connection. This fact simplifies some of the calculations below. 

Now, let {a;' 1 } and {x^ = x M — £^(x)} be two sets of coordinate systems, where ^(x) is an arbitrary — but 
infinitesimal, i.e., in this article, of first order — vector field. Then the components 4 Q '"' 3 fI ...„(i) of the tensor A with 
respect to the new coordinates x M can be related to the components of the tensor A a ' " s ^... v {x), defined with respect 
to the old coordinates with the help of the Lie derivative. Up to and including terms containing first order 

derivatives one has 

A a - 3 t ,..„{x) = A a ~%... v (x) + {CtA) a - p „...„(*) + ■ • • . (10) 

For a derivation of this equation, see Weinberg [14], Chap. 10, Sec. 9. 

Note that x in the left-hand side corresponds to a point, P say, of space-time with coordinates x^ in the coordinate 
frame {a;}, while in the right-hand side x corresponds to another point, Q say, with exactly the same coordinates x^ , 
but now with respect to the coordinate frame {x}. Thus, equation (10) is an expression that relates one tensor field 
A at two different points of space-time, points that are related via the relation (3). 

The following observation is crucial. Because of the general covariance of the Einstein equations, they arc invariant 
under general coordinate transformations x — > x and, in particular, under coordinate transformations given by (3). 
Hence, if some tensorial quantity A(x) of rank n (n = 0, 1, . . .) satisfies the Einstein equations with as source term 
the energy-momentum tensor T, the quantity A(x) = A{x) + C^A(x) satisfies the Einstein equations with source 
term T{x) = T(x) + C^T(x), for a universe with precisely the same physical content. Because of the linearity of the 
linearized Einstein equations, a linear combination of any two solutions is also a solution. In particular, C^A, being the 
difference of A and A, is a solution of the linearized Einstein equations with source term C(T . In first order, C^A(x) 
may be replaced by C^A m (t), where A m (t) is the solution for A(t) of the zero order Einstein equations. The freedom 
to add a term of the form £^A (0) (t), with (/i = 0, 1, 2, 3) four arbitrary functions of first order, to any solution of 
the Einstein equations of the first order, is the reason that none of the first order solutions is uniquely defined, and, 
hence, does not correspond in a unique way to a measurable property of the universe. This is the notorious gauge 
problem referred to in the introduction of this article. The additional terms C^A (0) (t) are called 'gauge modes'. 

Combining (9) and (10) we have 

A" ^ ^--.y (x) = A" ^ ^■■■ L i(x) + A" ^ '^•■■iy; T £ T 

AT---0 fO AOf-T 

)Li"'l/S \T ' ' ' -f 1 /^'"I^S \T 

+ A'*- 3 T ...„f ;lt + ■■■ + A a - : \... T ?, V . (11) 

We now apply the equation (11) to the case that A is a scalar cr, a four-vector and a tensor A^ v , respectively, 

a(x) = a{x) +C{x)d T a(x), (12a) 
V» = v" + V-rC - V T ^; T , (12b) 
Api, = A^v + A^vtSJ + A TV (J ' + A^ lT ^ T . v . (12c) 

For the metric tensor, we find in particular, from Eq. (12c), 

9nv = 9^u + in;u + (13) 



where we have used that the covariant derivative of the metric vanishes. 
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Our construction of gauge-invariant perturbations totally rest upon these equations for hatted quantities. In case 
<j(x) is some scalar quantity obeying the Einstein equations, <r(x) can be split up in the usual way in a zero order and 
a first order part: 

a(x) := <r m (t) +ct ( i)(x), (14) 
where cr m (t) is some background quantity, and hence, not dependent on the spatial coordinates. Then (12a) becomes 



&(x) =a m (t)+a (1) (x)+^(x)do(T m (t)+^(x)d ll a {1) (x). (15) 

The last term, being a product of the first order quantity ^(x) and the first order quantity d^a {l)1 will be neglected. 
We thus find 

a(x) = <r m (t) + a {1) (x), (16) 

with 

a {1) (x) := <T (1) (x) + tp(x)d a w (t), (17) 
where we used (5). Similarly, we find from (12b) and (13) 

v&(x) = v&+v& ;T e-vz,e-,r, as) 

and 

9wnv(x) = g m » v {x) + + (19) 

The latter two equations will be used later. 

We are now in a position that we can conclude the proof of the statement that ef^ and nf^ are gauge-invariant. 
To that end, we now write down the equation (17) once again, for another arbitrary scalar quantity u(x) obeying the 
Einstein equations. We then find the analogue of Eq. (17) 

uj(i)(x) = uj (1) (x) + il>(x)dauj m (t). (20) 

The left-hand sides of (17) and (20) give the value of the perturbation at the point with coordinates x with respect to 
the old coordinate system {x}] the right-hand sides of (17) and (20) contains quantities with the same values of the 
coordinates, x, but now with respect to the new coordinate system {x}. Eliminating the function tp(x) from Eqs. (17) 
and (20) yields 

- / x do<7 (0) (t) „ . d <7 (0 )(t) 

aM ~ d^itf M = aM - d^rf^ x) - (21) 

In other words, the particular linear combination occurring in the right-hand side of (21) of any two quantities lu and a 
transforming as scalars under an arbitrary space-time transformation — of first order, cf. Eq. (3) — is gauge- invariant, 
and, hence, a possible candidate for a physical quantity. 

The equation (21) is the key equation of this article as far as the scalar quantities e (1) and are concerned. It tells 
us how to combine the scalar quantities occurring in the linearized Einstein equations in such a way that they become 
gauge independent. The equation (21) can be used to immediately derive the expressions (2) for the gauge-invariant 
energy and particle number densities. 

In fact, let U fl (x) be the four- velocity of the cosmological fluid. In the theory of cosmological perturbations [see 
Eqs. (170)-(171) and (173)-(174)] there will turn out to be only three scalars, namely 

e(x) = c- 2 T^(x)U^x)U v (x), (22a) 

n(x) = c- 2 N^{x)U^{x), (22b) 

9{x) = c^U^x), (22c) 

where 



:= nW, 



(23) 
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is the cosmological particle current four- vector normalized according to U^U^ = c 2 . These scalars are split up 
according to 

e(x) = £((,)(*) +£(!)(#), (24a) 

n(x) = n m (t) +n m (x), (24b) 

8{x) = 9 m (t) + 6 w (x), (24c) 

where the background quantities e m (t), n m (t) and (o) (t) are solutions of the unperturbed Einstein equations. They 
depend on the time coordinate t only. The relation (21) inspires us to consider the gauge-invariant combinations 

e*(s) := e m (x)-^A Ull) (x), (25a) 

nUx) := n w (x)-^p^- Uw (x), (25b) 

:= (1) (x)-M^|c (1) (x). (25c) 

The question remains what to choose for u> in these three cases. In principle, for uj we could choose any of the following 
three scalar functions available in the theory, i.e., we could choose e. n or 6. As follows from (25a) and (25b), the 
choices u> = e and lu — n would lead to efl } (x) = or nfl } (x) = 0, respectively. This would mean that perturbations in 
the energy or the particle number density would show up only in the second order. Apparently, these choices would 
not be suitable since the perturbations of precisely these quantities are the object of study. We therefore are left with 
the choice 

u = 6, (26) 

which implies the equations (2a) and (2b) for the energy and particle number density perturbations, as was to be 
shown, and, moreover, 

0fl } = 0. (27) 

The latter equation implies that perturbations in the divergence of the cosmological velocity field 9(x) will only show 
up in second order. This constitutes no problem, since we are not interested in this quantity. Inserting (26) into (25a) 
and (25b) we obtain the expressions (2). Hence, it now has been shown that efl } and n^j, are indeed invariant under 
the infinitesimal coordinate transformation (3), i.e., that they are gauge-invariant. 

The linearized Einstein equations we are about to study contain the gauge function tp(x). If we eliminate s (1) (x) 
and n m {x) from these equations in favor of efi)(a;) and nfl } (x), with the help of the expressions (2), we obtain 
equations, which do not depend on ip(x) at all: all terms with ip(x) cancel automatically. This could be expected, 
because the equations should have solutions for ef^(x) and nfl^x) which are, by construction, independent of the 
diffcomorphism (3). Stated differently, the fact that ip(x) disappears from the linearized perturbation theory is in 
line with the earlier observation that (2) are gauge independent combinations for arbitrary diffcomorphisms (3). In 
Sec. VIII we will show in detail how this will happen. 



H. Non-relativistic limit of the gauge-invariant combinations 

The following question may arise. It is clear now that £ (1) (a;) and n {l) (x) are not gauge-invariant, but that the 
combinations sf\(x) and nfl\(x) are so. Hence, s {1) (x) and n {1) (x) cannot be physical perturbations, while efl^x) and 

nf 1 - ) (x) could be physical. But how can one be sure that the particular combinations (2) are indeed the real physical 
perturbations? The fact that they are gauge-invariant is a necessary, but not a sufficient reason. Since, moreover, 
any linear combination of gauge-invariant quantities is a new gauge- invariant quantity, it is not clear at all, yet, that 
the particular combinations (2) are the right ones. 

This issue can be settled by considering the non-relativistic limit, i.e., the limit of low spatial velocities with respect 
to the speed of light. In this particular case the linearized Einstein equations (202)-(203) imply, see Eqs. (228) and 
(229), 

VVfc) = 47rG £fl){t *' x) , (28) 
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where t p indicates the present time, c is the speed of light. G Newton's gravitational constant and V 2 is the usual 
Laplace operator. This is the well-known Poisson equation of the Newtonian theory of gravity. Hence, in the non- 
relativistic regime, the mathematical combination ef^ (2a), divided by c 2 , is to be interpreted as the normal, or 
'right' physical perturbation of the mass density g {1) , implying that e^ 1 , may be interpreted as the 'right' or physical 
perturbation in the energy density indeed. Since the non-relativistic limit of another Einstein equation, Eq. (220b), 
will turn out to imply nfl } (t p ,x) = efl } (t p ,x)/(mc 2 ), Eq. (226b), we see that also nf'j (2b) is the right physical 
perturbation. This answers the question raised at the beginning of this section with respect to ef^ and nfl } : these 
quantities are, indeed, the physical perturbations we were looking for. 



II. THE EINSTEIN EQUATIONS REWRITTEN IN SYNCHRONOUS COORDINATES 

As most authors do, we choose a synchronous system of reference. A synchronous system of reference is a system 
in which the line element for the metric has the form: 

ds 2 = c 2 d< 2 - 3ij (t, x)dx i dx j . (29) 

The name synchronous stems from the fact that surfaces with t = constant arc surfaces of simultaneity for observers at 
rest with respect to the synchronous coordinates, i.e., observers for which the three coordinates x % (i — 1, 2, 3) remain 
constant. A synchronous system can be used for an arbitrary space-time manifold, not necessarily a homogeneous or 
homogeneous and isotropic one. In a synchronous system, the coordinate t measures the proper time along lines of 
constant x 1 . From (29) we can read off that (x° = ct): 

goo(t,x) = 1, goi(t,x) = 0. (30) 

From the form of the line element in four-space, Eq. (29), it follows that minus gij(t, x), (i = 1, 2, 3), is the metric of 
a three-dimensional subspaces with constant t. Because of (30), knowing the three-geometry in all hypcrsurfaces, is 
equivalent to knowing the geometry of space-time. The following abbreviations will prove useful when we rewrite the 
Einstein equations with respect to synchronous coordinates: 

H ij :=-y ij , x l j:=g lk x k3 , H ij := +\g ij , (31) 

where a dot denotes differentiation with respect to x° = ct. From Eqs. (30)-(31) it follows that the connection 
coefficients of (four-dimensional) space-time 

r fiv = 29 P (dpfj.,^ + 9pu,t± ~ gfiv.p) ? (32) 

in synchronous coordinates are given by 

r°oo = r'oo = r°io = r°o; = o, (33a) 

r°ij = Hij, T\j = T l j = —>c l j, (33b) 
r fc y = \g kl {g H ,j + fa,., - gij,i) . (33c) 

From Eq. (33c) it follows that the T k ij are also the connection coefficients of (three-dimensional) subspaces of constant 
time. 

The Ricci tensor i? M „ := i? A M Ai/ is, in terms of the connection coefficients, given by 

R^lu r \_lv.\ r iix.v -i- r \ a r V(J . (34) 

Upon substituting Eqs. (33) into Eq. (34) one finds for the components of the Ricci tensor 

(35a) 
(35b) 
(35c) 

where the vertical bar in Eq. (35b) denotes covariant differentiation with respect to the metric gtj of a three-dimensional 
subspace: 

k 1 m := x* i|fc + V lk ^ 3 - T l jk x\. (36) 



R()i 


- k k k - 

= x k k\i 






Rij 


= Hij — 




t~ Rij j 
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The quantities 3 i?y in Eq. (35c) are found to be given by 

3r> T~*k j^k i j^l j^k j^l r^k /'Q7\ 

tiij = I — 1 ik.j ~r i yi fcZ — 1 ifci jl- (61) 

Hence, 3 Rij is the Ricci tensor of the three-dimensional subspaces of constant time. For the components i? M „ = g^ T R TV 
of the Ricci tensor (35), we get 

R° = x k k -x\x k u (38a) 

R\ = x k k]l ^x k l]k , (38b) 

R l j = x l 3 - x l jx k k + 3 R l j, (38c) 

where we have used Eqs. (30)-(31). 
The Einstein equations read 

where G^ v , the Einstein tensor, is given by 

G^ = R? v - \R a a g^. (40) 

In (39) A is a positive constant, the well-known cosmological constant. The constant k is given by 

^G . . 

k:=— 5-, 41 

with G Newton's gravitational constant and c the speed of light. In view of the Bianchi identities one has 

G^-u = 0, (42) 

hence, since g^ v M — 0, the source term of the Einstein equations must fulfill the properties 

T^ v . v = 0. (43) 

These equations are the energy-momentum conservation laws. An alternative way to write the Einstein equations (39) 
reads 

R% = K {T» V - \5^T a a ) - A£V (44) 

Upon substituting the components (38) into the Einstein equations (44), and eliminating the time derivative of x k k 
from the i?°o-cquation with the help of the i? l j-equations, the Einstein equations can be cast in the form 

{x k k f - 3 R- x k t x l k = 2(kT° + A), (45a) 
**kH " **i|k = nT t , (45b) 



x l 3 x k k + 3 R l 3 = k(T 1 3 - l^-T^p) - AS l j, (45c 



where 

3 R:=g ij3 R ij = 3 R k k , (46) 

is the curvature scalar of the three-dimensional subspaces of constant time. The (differential) equations (45c) are the 
so-called dynamical Einstein equations: they define the evolution (of the time derivative) of the (spatial part of the) 
metric. The (algebraic) equations (45a) and (45b) are constraint equations: they relate the initial conditions, and, 
once these are satisfied at one time, they are satisfied automatically at all times. 

The right-hand side of Eqs. (45) contain the components of the energy momentum tensor T Mt/ , which, for a perfect 
fluid, are given by 

T» v = (e+p)u»u„-p5» u , (47) 

where u^fx) = c~ 1 U ll (x) is the hydrodynamic fluid four-velocity normalized to unity: (u^u^ = 1), e(x) the energy 
density and p(x) the pressure at a point x in space-time. In this expression we neglect terms containing the shear and 
volume viscosity, and other terms related to irreversible processes. The equation of state for the pressure 

p = p(n,e), (48) 
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where n(x) is the particle number density at a point x in space-time, is supposed to be a given function of n and s 
(see also Appendix A for equations of state in alternative forms). 

As stated above already, the Einstein equations (45a) and (45b) are constraint equations to the Einstein equations 
(45c) only: they tell us what relations should exist between the initial values of the various unknown functions, in 
order that the Einstein equations be solvable. In the following, we shall suppose that these conditions are satisfied. 
Thus we are left with the nine equations (45c), of which, because of the symmetry of gij, only six are independent. 
These six equations, together with the four equations (43) constitute a set of ten equations for the eleven (6 + 3 + 1 + 1) 
independent quantities g^, u\ e and n. The eleventh equation needed to close the system of equations is the particle 
number conservation law iV+p = 0, i.e., 

(mi").,, = 0, (49) 

[see Eq. (23)] where a semicolon denotes covariant differentiation with respect to the metric tensor g^ v . This equation 
can be rewritten in terms of the fluid expansion scalar defined by Eq. (22c). Using Eqs. (33), we can rewrite the 
four-divergence (22c) in the form 

9 = u°- >c k k u" + S, (50) 

where the three-divergence $ is given by 

0:=u fc ]fc . (51) 

Using now Eqs. (22c), (33), (36) and (50), the four energy- momentum conservation laws (43) and the particle number 
conservation law (49) can be rewritten as 

f 00 + T ok lk +x k iT l k - x k k T 00 = 0, (52a) 

f l ° + T k \ k -2x l k T k0 - x k k T i0 = 0, (52b) 

and 

nu° + n tk u k + nO = 0, (53) 

respectively. Since T° l is a vector and T y is a tensor with respect to coordinate transformations in a subspace of 
constant time, and, hence, are tensorial quantities in this three-dimensional subspace, we could use in (52) a bar to 
denote covariant differentiation with respect to the metric gij(t,x) of such a subspace of constant time t. 

The Einstein equations (45) and conservation laws (52) and (53) describe a universe filled with a perfect fluid and 
with a positive cosmological constant. The fluid pressure p is described by an equation of state of the form (48): in 
this stage we only need that it is some function of the particle number density n and the energy density e. 

It is now time to actually derive the zero and first order Einstein equations. To that end, we expand all quantities 
in the form of series, and derive, recursively, equations for the successive terms of these series. Furthermore, we will 
now limit the discussion to a particular class of universes, namely the collection of universes that, apart from a small, 
local perturbation in space-time, are homogeneous and isotropic, the so-called Fricdmann-Lcmaitre-Robertson- Walker 
(flrw) universes. 



III. ZERO AND FIRST ORDER EQUATIONS FOR THE FLRW UNIVERSE 

Consider space-times which are a continuous collection (foliation) of three-dimensional, space-like slices of space- 
time, each of which is maximally symmetric. Expressed in synchronous coordinates, this statement means that we 
consider spaces with metric g^ v given by 

g 00 (t, x) = l, g 0i (t, x) = 0, (54a) 
gij(t, x) = -a 2 {t)g tJ {x), (54b) 

where gij(x) is the metric of such a three-dimensional maximally symmetric subspace. The minus sign in (54b) has 
been introduced in order to switch from the conventional four-dimensional space-time with signature (+,—,—,—) 
to the conventional three-dimensional spatial metric with signature (+,+,+). We write a 2 (t) rather than a{t) for 
the time dependent proportionality factor in (54b), the so-called scale factor or expansion factor of the universe, for 
reasons that will become clear later: the scale factor a(t) will turn out to be identifiable in some cases, as the 'radius of 
the universe', i.e., the slices of the foliation are three-dimensional surfaces in four-space, with radius a(t): see Eq. (77). 
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Essentially, however, we only suppose that the three-part of the metric can be factorized according to (54b), also in 
case the slices are no surfaces of hypcrsphcres. 

We now expand all quantities concerned in scries. We will distinguish the successive terms of a series by a subindex 
between brackets. 



£ 


— £(o) 4 


4t7 2 £ (2) H , 


(55a) 


n 


= n m ^ 


- V n d) 4 ?7 2 ^(2) H • 


(55b) 


P 


= P(o) 4 


- VPm + r l 2 P(2) H , 


(55c) 





= 0(o) 4 


■ rfi m 4 if9 {2) + • • • , 


(55d) 





= 0(„) ^ 


- rr& M 4 ?y 2 ^ (2) 4 ■ ■ ■ , 


(55e) 


''B 




+ ?7 3 i? (1) +ry 23 i?, 2 )4---- , 


(55f) 



where the subindex zero refers to quantities of the unperturbed, homogeneous and isotropic flrw universe. In order 
to derive the background and first order Einstein equations, we need ancillary quantities, which will also be expanded 
in series: 

m m = uf ) 4 r/itfi, 4- ifuf 2) H , (56a) 

9ij = 9(oyij + V9nnj 4 rfg^ij H , (56b) 

x%j = K {Q)ij 4 T}*C(i)ij + V 2>c (2)ij H , (56c) 

r"„ = r£,„ 4 rfr& v + tj 2 t£,„ + ■■■ , (56d) 

3 Rij = z R {a)lJ 4 n 3 R (1)ij + rj 2 % 2)ij + • ■ • , (56e) 
T ij = T\ 0) ij + V^(i)ij 4 r] 2 T^ 2)i j + • • • . (56f) 

In Eqs. (55) and (56) r\ (77 = 1) is a bookkeeping parameter, the function of which is to enable us in actual calculations 
to easily distinguish between the terms of different orders. 



A. Zero order quantities 



This section is concerned with the background or zero order quantities occurring in the Einstein equations. All 
results of this section are totally standard, and given here only to fix the notation unambiguously 

A tensor in a maximally symmetric space is called maximally form-invariant if its Lie-derivatives with respect to all 
Killing vectors of the maximally symmetric space vanish. Essentially, this means that the tensor is 'the same' in all 
directions and at all places, just as the metric of a maximally symmetric space is 'the same' in all directions and all 
places. It can be proved that the only maximally form-invariant scalar function which exists in a maximally symmetric 
space is a constant, that the only maximally form-invariant vector in a maximally symmetric space is the null vector, 
and that the only maximally form-invariant second rank tensor in a maximally symmetric space of dimension N > 3 
is proportional to the metric tensor (see Weinberg [14], Chap. 13). In our case, of a four-dimensional space-time which 
is a foliation of three-dimensional maximally symmetric spaces, this leads to the following, general, conclusions. A 
scalar can only be a function of time and the spatial part A % of a vector is the null vector. Furthermore, the spatial 
components F u of a tensor F^ 11 must be proportional to the metric g v and the components F ^ are identically zero. 

In particular, the background energy density e ia) (t,x), Eq. (55a), a space-time scalar, is independent of the coordi- 
nates x = (x , x , x ) of the maximally symmetric subspacc. Similarly, the background particle number density n (0) , 
Eq. (55b), the background pressure p m , Eq. (55c), and the fluid expansion scalar 6 m , Eq. (55d), all space-time scalars, 
are functions of time only. Thus we have, with respect to our particular coordinates, 

e (0) (t,x) = e (0) (i), n m (t,x) = n m (t), p (0) (t, x) = p (0) (t), 6 m {t,x) = 6 {0) (t). (57) 

The three-divergence $ (0) , Eq. (55e), and the curvature of the three-dimensional subspaces 3 i? (0 ), Eq. (55f), which are 
scalars only with respect to spatial coordinate transformations, are also functions of time only, i.e., 

i?(o)(* ) aj)=^(o)(*) ) % 0) (t,x)=% 0) {t). (58) 

Furthermore, the component of the four- vector u M , Eq. (56a), tangent to a maximally symmetric space is a three- 
vector in that maximally symmetric space. The only maximally form-invariant vector in a maximally symmetric space 
is the null vector. Hence, 



u\ a) = 0. 



(59) 
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Consequently, u M is proportional to S^q. Moreover, since w M is a unit vector, we have 

<„> = <*V (60) 

The zero order background metric tensor, occurring in Eq. (56b), has been supposed to be that of a maximally 
symmetric space, i.e., 

g mij (t,x) = -a 2 (t)g ij (x), (61) 

compare Eq. (54b). The time derivative of the three-part of the metric g^ij, K (a)iji Eq. (56c), may be expressed in 
the usual Hubble function H(i) := (da/ dt) / a(t) . We prefer to use a function H(t) = c _1 W(4), which we will call 
Hubble function also. Recalling that a dot denotes differentiation with respect to ct, we have 

H := -. (62) 
a 

Substituting the expansion (56b) into the definitions (31), we obtain 

x {0) ij = -Hg wi j, x* 0)j = -H5 t j , = -Hg^, (63) 

where we considered only terms up to the zero order in the bookkeeping parameter rj. 

Similarly, with Eqs. (50), (51), (55d), (55e), (60) and (63) we find for the background fluid expansion scalar, 8 {0) , 
and the three-divergence, $ (0) , 

m = 3H, d w = 0. (64) 

Using Eqs. (47), (55a), (55c), (56b), (60) and (61) we find for the components of the energy momentum tensor, 
Eq. (56d), 

^(0)0 = £ (0)i ^(0)0 = 0> ^To)j = ~P(°)S l j, (65) 

where the background pressure p m is given by the equation of state (48), which, for the background pressure, is 
defined by 

=p(n i0) ,e m ). (66) 

i.e., the subindex zero in p {0) refers to the zero order quantities it depends on; it is not a different function of its 
arguments. 

Finally, the Ricci tensor of the three-dimensional maximally symmetric subspaccs is proportional to the metric 
tensor of that subspace, i.e., 

3R (onj = Kgij, (67) 

(see Weinberg [14], Chap. 15, Sec. 1). The quantity K is time independent, as may be seen as follows. As follows 
from Eqs. (56e) and (56f), the background three-dimensional Ricci tensor, (37), is given by 

3p j^k pfc I pi j^k y\l pfc i ao\ 

■"-(ojij — 1 (o)ij,k 1 (o)ik,j ' 1 (o)ij L (o)kl 1 (a)ik L (o)jh \ {JO ) 

where the connection coefficients , are given by 



^(ojij — 2^(0) (5(o)K,j + 9(o)lj,i - 9(p)ij,l) , (69) 

where g 1 ^ and g@)ij depend on time. From Eq. (61), we have 

g g(t,x) = -J^g i i(x). (70) 
Hence, the connection coefficients r* 0) ^ are equal to the connection coefficients T k ij of the metric gif 

r f )ij = r fe y := \~g u (gu.j + gij,i - <%,/) ■ (7i) 
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Therefore, they do not depend on time. As a consequence, 3 R( )ij is time independent, implying that K is a constant. 
Maximally symmetric spaces may have curvature K that is positive, negative or zero. In the latter case the space is 
called flat. 

We now choose coordinates (r, 9, <f) such that the metric coefficients g m tj get the well-known Robertson- Walker 
form 

1 



1 r 2 

Or 2 sin 2 i 

where k = 0,±1. Comparing (61) and (72) we see that for this choice of coordinates we have 

9ij = diag { jT^ yy sin2 ■ ( 73 ) 

Substituting (72) into (68), combined with (71), we find 

3 R l0)ij = 2k gij . (74) 



We thus find, for the chosen rw coordinates 



From Eqs. (74) and (75) we have 



K = 2k. (75) 



Ik 



implying that the zero order curvature scalar 3 R( ) — 5?o) 3 R(anj is given by 

Pik 

%)(«) = (77) 

where we have used Eqs. (70) and (74). Note, that in view of our choice of the metric (+, — , — , — ), spaces of positive 
curvature k (such as spheres) have a negative curvature scalar 3 R (0 )- 

The results (61), (63), (65) and (74) are in agreement with the general observation that maximally form-invariant 
tensors of the second rank are proportional to the metric tensor of the space concerned. 

Thus we have found all background quantities: they are either space independent or proportional to g m ij(t,x), 
Eq. (61). The latter is proportional to (jij(x), the metric characteristic for spatial sections of constant time. 



B. First order quantities 



In this quite technical section we express all quantities occurring in the Einstein equations in terms of zero and first 
order quantities. The equations of state for the energy and pressure, e(n,T) and p(n,T), are not specified yet. 

Upon substituting the series (56a) into the normalization condition u^u^ = 1, one finds, equating equal powers of 
the bookkeeping parameter 77, 

<> = 0, (78) 
for the first order perturbation to the four-velocity. Writing the inverse of (56b) as 

ff W =<+wS + "-, (79) 
where is the inverse, (70), of g^M, (72), we find 

ffw = -9$)9m9wij, (80) 

and 

9wi = -9m9mli- (81) 
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It is convenient to introduce 

hij ■= ~g(i)ij, (82) 

so that 

W=g%, h i J = g^ ) h kj . (83) 
For the time derivative of the first order perturbations to the metric, Eq. (31), we get 



(i)y — 2 v> H mj ~ 2 i 1 K d) ~ 2 (^) 

The first order perturbation to the fluid expansion scalar 8, Eq. (50), can be found in the same way. Using (55d) 
and (60) one arrives at 

W = - \h\, (85) 
where we used Eqs. (78) and (84). The first order perturbation to the three-divergence Eq. (51), is 

^) (86) 

where we have used that 

(u fe | fc ) (1) = u*,| fc , (87) 

which is a consequence of 

r^fcu'eo) = 0. (88) 

The latter equality follows from Eq. (59). 

Upon substituting the series expansion (55a), and (55c)-(56b) into Eq. (47) and equating equal powers of r], one 
finds for the first order perturbation to the energy-momentum tensor 

J (l)0 — fc (i)> 

T (D0 = ( £ <o) +P(o)X), (89) 

where we have used Eqs. (60) and (78). The first order perturbation to the pressure is related to e {1) and n (1) by the 
first order perturbation to the equation of state (48). We have 

P(D =P n n m +p e e ( i), (90) 
where p n and p e are the partial derivatives of p(n, e) with respect to n and e, 



Since we consider only first order quantities, the partial derivatives are functions of the background quantities only, 
i.e., 

Pn = Pn[n i0) ,e m ), p e =p £ (n l0) ,s w ). (92) 

Using Eq. (33c) and the expansions (56b) and (56f), we find for the first order perturbations of the connection 
coefficients 

rfi)^ = —g^gwimT^yij + §<? (0) (g^uj + gwij,i — gayijj) ■ (93) 

The first order perturbation r*^-, Eq. (93), occurring in the non-tensor T k ij, Eq. (56f), happen to be expressible as 
a tensor. Indeed, using Eq. (82), one can rewrite Eq. (93) in the form 

r a)y = -3$) (&K| j + h W ~ h w)- (94) 
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Using the expansion for 3 Rij, (56e), and 1,^, (56f), one finds for the first order perturbation to the Ricci tensor (37) 

3d pfc pA: , pi pfc , pi pfc pi pA pi pA 

- 1 - (i)ij.k x (i)ik,j ' x (o)ij (i)fei ~r A (o)fei x (o)jfe i (i)»k A (o)j'i) \ aiJ ) 

which can be rewritten in the compact form 

3 R{i)ij = ^(i)ij\k ~ ^(i)ik\j- (96) 

By substituting Eq. (94) into Eq. (96), one can express the first order perturbation to the Ricci tensor of the three- 
dimensional subspace in terms of the perturbation to the metric and its covariant derivatives 

= ~\9{ )Q l ii\j\k + h-ij\i\k — hij\i\k — hjfc|i|j-)- (97) 

The perturbation ^R^j is given by 

*B* wj := (ff* 3 i? PJ )u) = 9% %»pi + | %o ) h i i , (98) 
where we have used Eqs. (61), (74), (77) and (83). Upon substituting Eq. (97) into Eq. (98) we get 

= — 5ff(o)(^ p\j\k + h k 3\v\k ~ h k\ P \j) + 55 ( o)^i|fe|i + | 3 R(o)h l j. (99) 

Taking i = j in Eq. (99) and summing over the repeated index, we find for the first order perturbation to the curvature 
scalar of the three-dimensional spaces 

3 R {1) = gg(h k mj - h k imk ) + § 3 R m h k k . (100) 

We thus have expressed all quantities occurring in the relevant dynamical equations, i.e., the system of equations 
formed by the Einstein equations combined with the conservation laws, in terms of zero and first order quantities to 
be solved from these equations. In the Sees. Ill C and III D below we derive the background and first order evolution 
equations, respectively. To that end we substitute the series (55) and (56) into the Einstein equations (45) and 
conservation laws (52) and (53). By equating the powers of rj° , r) , . . . , we obtain the zero order, the first order and 
higher order dynamical equations, constraint equations and conservation laws. We will carry out this scheme for the 
zero and first order equations only. 



C. Zero order equations 

With the help of Sec. Ill A and the expansions (55) and (56) we now can find from the Einstein equations (45) and 
conservation laws (52) and (53) the zero order Einstein equations and the conservation laws. Furthermore, in view of 
the symmetry induced by the isotropy, it is possible to switch from the six quantities and the six quantities to 
the curvature 3 R (0) (t) and the Hubble function H(t) only. 



1. Einstein equations 



Upon substituting Eqs. (63) and (65) into the (0, 0)-component of the Einstein equations, Eq. (45a), one finds 

3H 2 - \ 3 R (0) = ks (0) + A. (101) 

The (0, i)-components of the Einstein equations, Eqs. (45b), are identically fulfilled, as follows from Eqs. (63) and (65). 
We thus are left with the six (i, j)-components of the Einstein equations, Eqs. (45c). In view of Eqs. (63), (65) and (76) 
we find that x\ s) \ = xf 0)2 = ^0)3; ^(oji = ^(012 = ^(o)3 an d 3 ^f >i = 3 Rfo)2 = 3 ^?o>3' whereas for i ^ j these quantities 
vanish. Hence, the six (z, j)-components reduce to one equation, 

H = -3H 2 + | 3 R (0) + i/s(e (0 , - p m ) + A. (102) 

In Eqs. (101) and (102) the background curvature 3 i? (0 ) is given by Eq. (77). It is, however, of convenience to determine 
this quantity from a differential equation. Eliminating a(t) from Eqs. (62) and (77) we obtain 

3 R ia) +2H 3 R {0) =0, (103) 

where the initial value 3 i? (0) (io) is given by 

3 R m (t ) = --oTTT, (104) 

a [t ) 

in accordance with Eq. (77). It should be emphasized that Eq. (103) is not an Einstein equation, since it is equivalent 
to Eq. (77). It will be used here as an ancillary relation. 
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2. Conservation laws 



Upon substituting Eqs. (63) and (65) into the 0-component of the conservation law, Eq. (52a), one finds 

e co) +3ff(£ (0) +p (0) ) =0, (105) 

which is the relativistic background continuity equation. The background momentum conservation laws (i.e., the 
background relativistic Euler equations) are identically satisfied, as follows by substituting (63) and (65) into the 
spatial components of the conservation laws, Eq. (52b). 

The background particle number density conservation law can be found by substituting (60) and (64) into Eq. (53). 
One gets 

h m + 3Hn m = 0. (106) 

This concludes the derivation of the background equations. 

What we have found is that for the Robertson- Walker metric (72) the Einstein equations (45) and conservation 
laws (52) and (53) reduce to the initial value condition (101) and four differential equations (102), (103), (105) 
and (106) for the four unknown functions H, 3 i?( ), £( ) and n (0) . If we use Eqs. (64) and (66), we can state that we 
have found equations for all unknown background quantities (55) in zero order approximation, namely (o) , i9 (0) , 3 i? (0) , 
e (0) and n (0) . In Sec. V we derive equations for the corresponding first order quantities 3 i? ( i), e ( i) and n (1) . 



D. First order equations 



In this section we derive the first order perturbation equations from the Einstein equations (45) and conservation 
laws (52) and (53). The procedure is, by now, completely standard. We use the series expansion in r\ for the 
various quantities occurring in the Einstein equations and conservation laws of energy-momentum, and we equate the 
coefficients linear in rj to obtain the 'linearized' or first order equations. 



1. Einstein equations 



Using the series expansions for 3 R, (55f), ?f*j, (56c), and T°q, (56d), in the (0, 0)-component of the constraint 
equation (45a), one finds 

2*?„, fc *ki - %r> - 2x k {0)l x\ 1)k = 2*2*0. (107) 

With the zero order equation (63), the abbreviation (84) and the expression for T® 1)0 , (89), we may rewrite this 
equation in the form 

Hh k k + ± 3 R {1) =-KE m . (108) 

Employing the series expansions for x l j, (56c), and T°i, (56d), we find for the (0, z)-componcnts of the constraint 
equations (45b) 

^W-4)iife = KT °)f ( 109 ) 

where we noted that 

(**j|k)w = *fi>,-|k> ( n °) 

which is a consequence of 

which, in turn, is a direct consequence of Eq. (63). From Eqs. (84) and (89) we find 

h k k]l - h k i\ k = 2k( £(0) + p i0) )u {l)l . (112) 

Finally, we consider the (i, j)-components of the Einstein equations (45c). Using the series expansions for >c l j, 
(56c), T l 3 , (56d), and 3 i?% , (56e), we find 
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With Eqs. (63), (84) and (89), we get 

h i j + 3Hti j + 5 i j Hh k k + 2 3 Rl )j = - K ^(e (1) (114) 

where 3 R\ 1} j is given by Eq. (99). 

Note that the first order equations (108) and (114) are independent of the cosmological constant A: the effect of 
the non-zero cosmological constant is accounted for by the zero order quantities [cf. Eqs. (101) and (102)]. 



2. Conservation laws 



We now consider the energy conservation law (52a). With the help of the scries expansions for (56c), and T^ v , 
(56d), one finds for the first order equation 

+ T^j \ k + >cf a) [T l (1)k + ^ )Z T ( ' 0)fe — j^ 0)k T^ — x (1)i .T ( °j = 0, (115) 
where we have used that for a three-vector T ok we have 

(T V) (1 ,=T^| fc , (116) 
see Eq. (87). Employing Eqs. (63), (65), (84), (86) and (89) we arrive at the first order energy conservation law 

i w +3H(e w +p w ) + (e t0 , + P(0) )« ) | fc - \h k k ) = 0, (117) 

Next, we consider the momentum conservation laws (52b). With the series expansions for (56c) and T Miy , 
(56d), we find for the first order momentum conservation law 

T'iO I (rpik \ r>~^ rpkO O-^ rpkO „_fc TniO „ Ji rpiQ f\ /i 1 n\ 

(i) + l J fcj(i) _ Z>C ( )k 1 (i) - ZK (l)k 1 (0) ~ X (0)fc J (l) ~ >C (l)k 1 (Q) ~ U - l ilS J 

Using that 

(T lk { kU = -9%P W \k, (H9) 

and Eqs. (63), (65), (84) and (89) we arrive at 

- 9%P w \h + 5Jr(e ( o) +P(o))< ) = 0, (120) 

where we have also used that the covariant derivative of g 1 ^ vanishes: <7*^i fc = 0. 

Finally, we consider the particle number density conservation law (53). With the expansions for n, (55b), 9, (55d), 
and it'*, (56a), it follows that the first order equation reads 

h m u° m + h m u° w + n^.fciifij + «(i),fcwf ) + n m G w + ?vAo) = 0. (121) 
With the help of Eqs. (60), (64) and (78) we find for the first order particle number conservation law 

n (1) + 3Hn w + n w (u k 1){k - \h k k ) = 0, (122) 
where we employed Eq. (85) to eliminate 



i a 



( e (0) +P(0))U(i) 



3. Summary 



In the preceding two Sees. Ill D 1 and III D 2 we have found the equations which, basically, describe the perturbations 
in a flrw universe, in first approximation. They are Eqs. (108), (112), (114), (117), (120) and (122). For convenience 
we repeat them here 



Hh k k 
h k 

n k. u. - 



f- 2 ^R(l) — ^£(1) i 

h k i\ k = 2k(£ ( o) +p {0) )u {1) 



h l 3 + 3Hh l + 8 l 3 Hh k k + 2 *R> , = -K6 l j(s (1) - p w ) 



£ {1) + 3H(e (1) 
1 d 



■P(o))(tt* 



)\k 



1 Lk 



h k k ) = 0, 



at 



( £ (0) 



5(o> ( i)|fe + 5tf(e (0) + p i0) )u l (1) = 0, 



n (1) + 3Hn (1) +ri ( o)Ki ) | fc - \h k k ) = 



(123a) 
(123b) 
(123c) 
(123d) 

(123e) 
(123f) 
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where 3 R\ t) j and 3 i? (1) are given by Eqs. (99) and (100), respectively. Hence, the equations (123) essentially are fifteen 
equations for the eleven (6 + 3 + 1 + 1) unknown functions h?j, u % m , e (1) and n {1) . The pressure p {0) is given by an 
equation of state (66), and the perturbation to the pressure, is given by Eq. (90). The system of equations is not 
overdetermined, however, since the four equations (123a) and (123b) are only conditions on the initial values. These 
initial value conditions are fulfilled for all times t automatically if they are satisfied at some (initial) time t = to. 



IV. CLASSIFICATION OF THE SOLUTIONS OF FIRST ORDER 



The set of equations (123), which are linear in their (eleven) unknown functions, can be split up into three sets 
of equations, which, together, are equivalent to the original set. We will refer to these sets by their usual names of 
scalar, vector and tensor perturbation equations. We will show that the vector and tensor perturbations do not, in 
first order, contribute to the physical perturbations ef^ and nfl } . As a consequence, we only need, for our problem, 
the set of equations which are related to the scalar perturbations. By considering only the scalar part of the full set 
of perturbation equations we are able to cast the perturbation equations into a set which is directly related to the 
physical perturbations efl } and n^ 1 ,. This is the subject of the next section, Sec. V. 

At the basis of the replacement of one set (123) by three sets of equations stands a theorem proved by Stewart [15, 16] , 
which states that a symmetric second rank tensor can be split up into three irreducible pieces, and that a vector can 
be split up into two irreducible pieces. Here, we will use this general theorem to obtain equations for the scalar 
irreducible parts of the tensors h l j and 3 R l {1) j and the vector namely hlj, 3 R\ 1) \\j and u (1) \\. 

For the perturbation to the metric, a symmetric second rank tensor, we have in particular 

tii=h\i+h±i+Ki, (124) 
where, according to the theorem of Stewart [15], the irreducible constituents hlj, h l j_j and h\j have the properties 

h\j = !(#%• + C |i | j )> (125a) 

hit = 0, (125b) 
hl k = 0, h*i\ h = 0, (125c) 

with <f>(t,x) and £(t, x) arbitrary functions. The contravariant derivative A' 1 is defined as g^Ay. The functions hij, 

h' l ± j and h\j correspond to scalar, vector and tensor perturbations, respectively. 

In the same way, the perturbation to the Ricci tensor can be decomposed into irreducible components, i.e., 

3 ^ j = 3 % ) \\ j + 3R U± j + 3R U* j - (126) 
The tensors ^jjiui 3 R(i)±j an d 3 ^u)*j nave the properties comparable to (125), i.e., 

3R U\\j = l^+^li), (127a) 
3 R^±k = 0, (127b) 
3 ^,* fc = 0, 3 ^ i|fc = 0, (127c) 

where j(t, x) and 7r(t, x) are two arbitrary functions. By now using Eq. (99) for each of the irreducible parts we find 



\j + S i j<t> lk \k~<; lkli \ j \k-( lk \ j li \k + { lk \k li \ j + C li \j lk \k + f %)(^ + C |1 |i)l, (128a) 



K W\\j ~ c 

3r Id±j = -jz9m(h k L P \j\k + + \9(p)h\j\k\i + \ 3R m^xj> (128b) 

3Rl (i)*j = ~l9lo)(h k P \j\k + ^*j|p|fe) + ^9( )Kj\k\l + | 3 R(o)Kj- (128c) 

Combining the expressions (127) and (128), we can derive relations, to be obeyed by 7, tt, h l ± j, and hlj. Firstly, from 
property (127a) and (128a) it follows that 

7 = W lk \k + l% 0) <P), (129a) 



7T 



= 5(^i,-C |fe!t i,| fe -C |fc i, |i | fe + C |A; | fc |j | J + C |i i, |fc | fc + | 3 R ( o ) C |l i,). (129b) 
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In a flat flrw universe, Eqs. (129) reduce to 



implying that, for a flat FLRW universe, 



7 = W k \k, (130a) 
71-1*1* = (130b) 



l = ^ k \k, (131) 



whereas there is no such restriction on the functions </> and ( in a flat flrw universe. 
Secondly, combining expressions (127b) and (128b) it follows that h\j has the property 

hf\ k \i = 0, (132) 

in addition to the property (125b). In Sec. IV B we show that this additional condition is needed to allow for the 
decomposition (141). 

Finally, combining Eqs. (127c) and (128c), we find that h\j must obey 

g k> '{h™ k\i\m\l + h™i\k\m\l ~ K?i\k\l\m) = 0, (133) 

in addition to (125c). The relations (133) are, however, fulfilled identically for flrw universes. This can eas- 
ily be shown. First, we recall the well-known relation that the difference of the covariant derivatives A 1 '" 1 k ...i\ p \ q 
and A l '" J k...i\ q \ p of an arbitrary tensor can be expressed in terms of the curvature and the tensor itself (Weinberg [14], 
Chap. 6, Sec. 5) 



A k---l\p\q 


A- 3 k---l\q\p 








3ds i 
•« (o)kpq ' 


■ + A i -\., 


3ps 
■s ''(o)ipg 


-A s -i k . 


*** ^(tyspq 


■-A*-%. 


i 3 R j 
■■I - LL (0)spqi 



(134) 

where R l ^j k i is the Riemann tensor for the spaces of constant time. At the right-hand side, there is a term with a 
plus sign for each lower index and a term with a minus sign for each upper index. 
We apply this identity taking for A the second rank tensor h\j to obtain 

urn _ j m urn 3ns _ i s 3pm /l ocl 

"•* k\i\m "* fc|m|i — "* s fl ([)) i lm n *k "(ojaim" \ 1JJ I 

Now note that h^ k \ m u vanishes in view of (125c). Next, we take the covariant derivative of (135) with respect to x l , 
and contract with g kl 

g kl h™ k \i\ m \i = g kl {h™ s 3 R s i0)kirn - K k 3 R™ )sim )\i- (136) 
Next, using the expression which one has for the Riemann tensor of a maximally symmetric three-space, 

3 % )bcd = k(6 a c g bd ~S a d g bc ), (137) 
(where k = 0,±1 is the curvature constant) we find 

g M h? kli]mll = 0, (138) 
i.e., the first term of (133) vanishes. The second and third term can similarly be expressed in the curvature 

urn i)7i i m 3ns i um 3ps _ us 3pm /i on\ 

"* i|fe|m|; z|fe|i|m — "■* s\k ■ n -(o)imi ' a * 'I s n (o)/;ml IL *i\k -"-(ojam/) \ 10J ) 

where the general property (134) has been used. Upon substituting the Riemann tensor (137) and contracting with g kl , 
we then arrive at 

g kl (h™ lWm{l -h? tmm )=0, (140) 

i.e., the second and third term of (133) together vanish. Hence, for flrw universes, Eq. (133) is identically fulfilled. 
Consequently, the decomposition (125c) imposes no additional condition on the irreducible part h\ j of the perturbation 
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The three-vector can be uniquely split up according to [15] 

= «(i)|| +u w ±, (141) 

where is the longitudinal part of u il} , with the properties 

VA(u {1 ,||) = 0, V-u (1) = V-u tl ,||, (142) 

and U(ijjl is the transverse part of u (1) , with the properties 

V-w (1) _l = 0, VAiip) = VAttp,!, (143) 

where the divergence of the vector tt (1) is defined by, Eq. (51), 

V ■«(!,:=«*,)* = (1) , (144) 

and the rotation of the vector tt (1) is defined by 

(V A u a) )i := ei jk u {1)j \ k = e l jk u (1)hk , (145) 

where e/ is the Levi-Civita tensor with ei 23 = +1. In Eq. (145) we could replace the covariant derivative by the 
ordinary partial derivative because of the symmetry of r*^. 

Having decomposed the tensors h l j, R-l^j an d in a scalar ||, a vector _L and a tensor part *, we can now 
decompose the set of equations (123) into three independent sets. The recipe is simple: all we have to do is to append 
a subindex || , _L or * to the relevant tensorial quantities in equations (123). This will be the subject of the Sees. IV A, 
IV B and IV C below. 



A. Tensor perturbations 

We will show that tensor perturbations are not coupled to, i.e., do not give rise to, density perturbations. 
Upon substituting h l j = h\j and 3 -R( 1} j = ^m*,- into the perturbation equations (123) and using the proper- 
ties (125c) and (127c), we find from equations (123a), (123b) and (123d) 

£ (1) = 0, p (1) = 0, n (1) = 0, it (1) = 0, (146) 

where we have also used Eq. (90). With (146), Eqs. (123e) and (123f) are identically satisfied. The only surviving 
equation is (123c), which now reads 

Kj I -Ulh;, ■ ■>■■■!!■ ., li. (147) 
where 3 i?j 1} ^ is given by Eq. (128c). Using Eqs. (85), (86), (125c) and (146) it follows from Eqs. (2) that 

4=0, n& = 0, (148) 

so that tensor perturbations do not, in first order, contribute to physical energy density and particle number density 
perturbations. Hence, the Eqs. (147) do not play a role in this context, where we are interested in energy density and 
particle number density perturbations only. 

The equations (147) have a wave equation like form with an extra term. The extra term 3Hh\j in these equations is 
due to the expansion of the universe. Therefore, these tensor perturbations are sometimes called gravitational waves. 
The six components h*j satisfy the four equations (125c), leaving us with two independent functions hlj. They are 
related to linearly and circularly polarized waves. 

B. Vector perturbations 

We will show that, just like tensor perturbations, vector perturbations are not coupled to density perturbations. 
Upon replacing h l j by h\- and ^R\ 1} j by 3 J2* 1J _ LJ - in the perturbation equations (123), and using the expressions (125b) 
and (127b), we find from Eq. (123a) and the trace of Eq. (123c) 



£(D = 0, P(d = 0, n (1) = 0, 



(149) 
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where we have also used Eq. (90). 

Since h\- is traceless and raising the index with g^ 0) in Eq. (123b) we get 

h k f ]k + 2Hh k l\ k = 2/s(e {0) +p {0 )H), (150) 

where we have used Eqs. (31) and (63). We now calculate the covariant derivative of (150) with respect to x J , and 
use (132) to obtain 

V ■«(!)=(), (151) 

where we made use of the fact that the time derivative and the covariant derivative commute. With Eqs. (141)-(143) 

we see that only the transverse part of u (lu namely plays a role in vector perturbations. From (125b) and (149) 
it follows that the equations (123d) and (123f) are identically satisfied. The only surviving equations are (123b), 
(123c) and (123e), which now read 

h±i\k = -2n(s m +p m )u m j_i, (152a) 

h\ 3 + ZHh^j + 2 3 Rl 1)±J = 0, (152b) 

' r (e ( o) +P m )u\ 1)± ] + 5H{e (a) +p m )u\ 1)± = 0, (152c) 



c 



dt 



where ^^j^ is given by Eq. (128b). 

Using Eqs. (85), (86), (149) and (151) we get from Eqs. (2) 

4=0, n&=0, (153) 

implying that also vector perturbations do not, in first order, contribute to physical energy density and particle 
number density perturbations. Hence, the equations (152) do not play a role when we are interested in energy density 
and particle number density perturbations, as we are here. Vector perturbations are also called vortices. 

Since vector perturbations obey V • u (1) ± = 0, they have two degrees of freedom. As a consequence, the tensor h\- 
has also two degrees of freedom. These degrees of freedom are related to clockwise and counter-clockwise rotation of 
matter. 



C. Scalar perturbations 

Differentiation of Eqs. (123b) covariantly with respect to x J we obtain 

h\k\i\j - h\i\k\j = 2/t(e C0 ) +P(o)) u (i))!i|r (154) 
Interchanging i and j in this equation, and subtracting the resulting equation from (154) we get 

~~ = ~ 2k(£(o) + P(o))(W(i)||i|j — U(i)\\j\i), (155) 

where we have used that ft|fe|j|j = ^|fc|j|i- Using that V Ait (1) || = 0, we find from Eq. (125a) that the function £ must 
obey the equations 

C |fc |i|fcU-C |fe U|fcK=0. (156) 

These equations arc fulfilled identically in flrw universes. This can be seen as follows. We first rewrite these equations 
by interchanging the covariant derivatives in the form 

(C |fc |i|fcii - C |fc |ibife) - (C |fc ii|fc|i - C |fe ii|i|fe) + (C |fc |i|j - C |fe |j|i)|fc = o. (157) 

Next, we use Eq. (134) and substitute the Riemann tensor (137) into the resulting expression. Using that = Ch'li) 
we find that the left-hand sides of the Eqs. (157) vanish. As a consequence, the Eqs. (156) are identities. Therefore, 
the decomposition (125a) imposes no additional condition on the irreducible part h l »j of the perturbation h l j. 
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The evolution equations (123) for scalar perturbations read 



Hh\ k 



1 3p 

2 -"(Dll 



-KS, 



h\\k\i — hp\k — 2/c(e co) +p m )u w \\i, 

e (1) + 3i/(e (1) +p (1) ) + (e (0) + P m )[{u k m p\ k - |/i| fc ] = 0, 



1 d_r 

cut 

n m +3Hn {1) + n m [(u^) lk - \h\ k 



(£<o) +P(o) Kdii - ff(o>d)|fc + 5i?(e (0) +p (0) )W( 1) || = 0, 



(158a) 
(158b) 
(158c) 
(158d) 

(158e) 
(158f) 



where the perturbations to the metric and the Ricci tensor are given by Eqs. (125a) and (128a), respectively. In the 
tensorial and vectorial case we found e (1) = and = 0, implying that efl } — and nf[^ = 0, which made the tensorial 
and vectorial equations irrelevant for our purpose. Such a conclusion cannot be drawn from the equations (158). 
Perturbations with e (1) ^ and n (1) ^ are usually referred to as scalar perturbations. 

In Sec. V we rewrite the set of equations (158) in such a way that they determine the evolution of the quantities (2). 
Since the perturbation equations contain only the components hp, it follows that relativistic energy density and 
particle number density perturbations are characterized by two potentials, (f> and £. 



D. Summary 

In the foregoing three Sees. IVA-IVC we have written down, using the decompositions (124), (126) and (141), the 
perturbation equations for gravitational waves, Eqs. (147), vortex perturbations, Eqs. (152) and scalar perturbations, 
Eqs. (158). For gravitational waves and vortices it was shown that ef^ = and nf[ } = 0. Therefore, we consider from 
now on only the Eqs. (158) for the scalar perturbations. 

An important consequence of the theorem of Stewart, Eqs. (124)-(125) and (126)-(127), is that we need only two 
potentials <j> and C, (instead of six potentials h?j) to describe the evolution of the scalar perturbations. In Sec. IX 
we will show that, in a flat flrw universe, the potential £ automatically disappears from the perturbed Einstein 
equations and that </>(£, x) can be related to the well-known, time-independent potential y>(x) encountered in the 
Poisson equation V 2 ip(x) = A-kGq {1) (x), where g m {x) = £^)(i P , x)/c 2 is the mass density at x. In fact, we will show 
that in the non-relativistic limit we have f{x) = <j>{x) / a 2 (t p ) , where a(t) is the scale factor occurring in (54b). In 
these equations t p stands for the present time. 

In the definitions (2) of the physical perturbations ef^ and nf^ the background quantities e (0) , n (0) and 6 {0) = 3H 
occur. The evolution of these quantities is governed by the background Einstein equations given in Sec. Ill C. Of the 
first order quantities e (1) , n (1) and 9 {1) in (2), only the first two do explicitly show up in the first order equations (158). 
In Sec. V we rewrite this set of equations in such a way that also the first order quantity 8 (1) does explicitly occur 
in the set of perturbation equations. The result is that the evolution of the physical quantities (2) is completely 
determined by the background- and first order equations. 



V. SCALAR FIRST ORDER EQUATIONS 

In this quite technical section we rewrite the scalar perturbation equations (158) in terms of quantities 9 {1) , 3 -R ( i)||, 
£(i) and which are suitable to describe exclusively the scalar perturbations. The result is that the first order 
quantities occurring in the definitions (2) do explicitly occur in the set of equations. 
Eliminating the quantity h^ k from Eq. (158a) with the help of Eq. (85) yields 

2!/(# (1) -tf (1) )-i 3 i? (1) || =«e (1) . (159) 

Thus the (0, 0)-component of the constraint equations becomes an algebraic equation which relates the first order 
quantities 6 {1) , 3 i? ( i)|| and e {1) . 

It now takes some steps to rewrite the three constraint equations (158b) in a suitable form. Firstly, multiplying 
both sides by g^ 0) and taking the covariant divergence with respect to the index j we find 

ff(0)(^j|Jfe|t|j ~ ^fi|A|j) = 2 «(S(o) +P(o)Ki), (160) 
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where we have also used Eq. (86). The left-hand side will turn up as a part of the first order derivative of the curvature 
3 i? (1) ||. In fact, differentiating Eq. (100) with respect to ct and recalling the fact that the connection coefficients ■, 
Eqs. (71), are independent of time, one gets 

3 7? (1) || = -2H 3 R (1)[l +g% (h\ m - h\ mo ) + \*R {a) h\ k , (161) 

where we have used Eqs. (31), (63), (103) and 

9(li h p\j\k = g\l)h\i\k\j, (162) 

which is a consequence of g l ^\ k = and the symmetry of h 1 ^ . Next, combining Eqs. (160) and (161), and, finally, 

■ 



eliminating /ifffc with the help of Eq. (85), one arrives at 



3 i? (1) H + 2H 3 i? (1) |, - 2/s(e (0) +pMd + f 3 i?<o)(#<i> - W ) = 0. (163) 

In this way we managed to recast the three (0, i)-components of the constraint equations in the form of one ordinary 
differential equation for the local perturbation, || , to the spatial curvature. 

We now consider the dynamical equations (158c). Taking the trace of these equations and using Eq. (85) to eliminate 
the quantity h\ik> we arrive at 

6 {1) - + 6H(9 W - tl) ) - 3 i? (1) |, - §«(e (1) - p (1) ). (164) 

Thus, for scalar perturbations, the three dynamical Einstein equations (158c) with i = j reduce to one ordinary 
differential equation for the difference 9 {l) — For i ^ j the dynamical Einstein equations are not coupled to scalar 
perturbations. As a consequence, they need not be considered. 

Taking the covariant derivative of Eq. (158e) with respect to the metric g^ij and using Eq. (86), we get 

^ (e<o) +P(o,)i? ( i) - g\o)Pw\k\i + 5H(e (0) + p (Q) )tf (1) = 0, (165) 

where we have used that the operations of taking the time derivative and the covariant derivative commute, since the 
connection coefficients r* 0)jj -, (71), are independent of time. With Eq. (105), we can rewrite Eq. (165) in the form 

# +H (2- 3^-) <? w + ] V \ w = 0, (166) 

where V 2 is the generalized Laplace operator which, for an arbitrary function f(t, x) and with respect to an arbitrary 
three-dimensional metric g tJ (x) 7 is defined by 

V 2 / := fl* J 7|i|r (167) 

Thus, the three first order momentum conservation laws (158e) reduce to one ordinary differential equation for the 
divergence $ ( i). 

Finally, we consider the conservation laws (158d) and (158f). Eliminating the quantity h^k from these equations 
with the help of Eq. (85), we get 

i m + 3H(e m +Pm) + (e w +Pm)0 w = 0, (168) 

and 

+ 3Hn (1) +n (0) 9 (1) = 0. (169) 

The algebraic equation (159) and the five ordinary differential equations (163), (164), (166), (168) and (169) is a 
system of six equations for the five quantities m , 3 i? ( D, s {1) and n (1) , respectively. This system is, however, 
not overdetermined, since the constraint equation (159) is only an initial value condition. It can easily be shown by 
differentiation of Eq. (159) with respect to time and eliminating the time derivatives with the help of Eqs. (101), 
(102), (163), (164) and (168) that the general solution of the system (159), (163), (166), (168) and (169) is also a 
solution of Eq. (164). Therefore, we do not consider Eq. (164) anymore and we take the algebraic equation (159) 
and the four ordinary differential equations (163), (166), (168) and (169) for the five quantities 6* (1) , 3 R W , £ {1) 
and n (1) , respectively, as the basis of our perturbation theory. For scalar perturbations this system of equations is 
equivalent to the full set (123) of linear Einstein equations and conservation laws. 
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VI. SUMMARY QUANTITIES AND EQUATIONS 
A. Zero order equations 

The Einstein equations and conservation laws for the background of an flrw universe are given by (102), (105), 
(106) and (103): 

(170a) 
(170b) 
(170c) 
(170d) 



H = 


~£ 3 ^(0) — 5 K£ (0)(1 + W )l 


£ (0) = 


-3He (0) (l + w), 


n (0) = 




-fl(O) — 


—2H 3 R (0) , 



and the constraint equation (101) 



3H 2 = 1 3 i? (0) + Ke m + A. (171) 



In order to arrive at Eq. (170a) we eliminated the term 3H 2 from Eq. (102) with the help of Eq. (101). Furthermore, 
we introduced the abbreviation 

«,(*) := P f^-. (172) 

E (0) W 

The set (170) consists of four differential equations with respect to time for the four unknown quantities e (0) (t), n {0) (t), 
(o) (t) = 3i/(t), and 3 R (a) (t). Recall # (0) = 0, Eq. (64). The pressure p^it) is related to the energy density £^{t) and 
the particle number density n (0) (t) via the equation of state (66). The algebraic equation (171) is a constraint on the 
initial values. 



B. First order equations 

The first order equations describing density perturbations are given by the set of four differential equations (168), 
(169), (166) and (163) 

e (1) +3H(e m +p w ) + e (0) (l + w)6 (1) = 0, (173a) 
n (1) + 3Hn w + n w 6 w = 0, (173b) 

^ (1) +g(2-3/3 2 K 1)+ 1 V ^ (1) ^0, (173c) 

3 # (1) |l + 2H 3 i? (1) |, - 2/se (0) (l + w)0 w + 1 3 i? ( o)(^ (1 ) - = 0, (173d) 
together with one constraint equation, Eq. (159) 

2tf(# (1) -tf (1) )-i 3 i? (1) || =«e (1)I (174) 

for the five unknown functions £ ( i), and # ( i), respectively. These are the first order perturbations to 

the background quantities e (0) , n (0) , $ (0) = 0, 3 i? (0 ), and 6 m = 3H. The first order perturbation to the pressure is 
given by the perturbed equation of state (90). 

As has been shown in Sec. V, the equations (173)-(174) comprise the conservation laws and the constraint equations. 
As a consequence, the general solution of these equations is also a solution of the dynamical Einstein equation (164). 
Therefore, we may, in the study of scalar perturbations, replace the full set of perturbation equations (158) by 
the set (173)-(174). Note that only three scalars e, n and 8 :— play a role in a density perturbation theory. 
The only non-trivial gauge-invariant combinations which can be constructed from these scalars and their first order 
perturbations are the combinations (2). These combinations do always exist, since, in a non-static universe, we have 
0(o) (t) 7^ for all times t, as follows from Eq. (170a). 

The operator V 2 , occurring in Eq. (173c), is the generalized Laplace operator defined by Eq. (167). The quantity (3(t) 
occurring in Eq. (173c) is defined by 

V £co)(i) 
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Eliminating the time derivatives of e (0) and n {0) from Eqs. (2) with the help of the background equations (170), we 
arrive at 



6He m (l + w) 
3 R m +3ne (0) (l + w)' 
i = _ 6Hn {0) 
(1) (1> 3 i? (0 ) +3Ke m {l + w) 



(176a) 
(176b) 



where we have used that m = 3i?, Eq. (64). We have achieved now that the set of perturbation equations (173) 
and (174), together with the background equations (170) and (171) determine the evolution of the physical quan- 
tities (2), or, equivalently, (176). In principle, we are ready. However, the solution of the set of equations (173) 
and (174) is gauge dependent (see Appendix B). Indeed, one may easily check, by a direct calculation, that these 
equations are invariant under the transformation 

(177a) 
(177b) 

(177c) 

2 l RwiA , (177d) 







£(i) + ty£(o)i 


n m - 


-> n (1 ) = 


n w + tph i0) , 




- tf« = 


a 2 ' 




-> 3 R(d\\ 


= 3 -R (1) ||+4J? 




-> 9 {1) = 


+'#(0), 



w 13 



(177c) 

where V 1 is time independent in synchronous coordinates, see Eq. (B2). By switching from the variables e (1) , n (1) , t9 (1) , 
3 i? (1) H and 9 {1) to the variables efl } and nf^, we will arrive at a set of equations for efl^ and nfl } with a unique, i.e., 
gauge-invariant solution. This will be the subject of Sec. VIII. First, however, we derive some auxiliary equations 
related to the entropy, pressure and temperature. 



VII. ENTROPY, PRESSURE AND TEMPERATURE 
A. Gauge-invariant entropy perturbations 



The second law of thermodynamics reads 

dE = TdS-pdV + (idN, (178) 

where E, S and N are the energy, the entropy and the number of particles of a system with volume V, and where /z, 
the thermodynamic — or chemical — potential, is the energy needed to add one particle to the system. In terms of 
the energy per particle E/N = e/n and the entropy per particle s = S/N the law (178) can be rewritten 

—N^j = T d(sN) -pd(^j+[idN, (179) 

where e and n are the energy and particle number densities. With the Euler relation /i = (e + p)n~ 1 — Ts, the second 
law can be cast in a form without \x and N. In fact, using the rule d(/<?) = fdg + gdf we immediately find from 
Eq. (179) 

Tds = d(^j +pd(^), (180) 

where p, and TV have canceled, indeed. The thermodynamic relation (180) is true for a system in thermodynamic 
equilibrium. For a non-equilibrium system that is 'not too far' from equilibrium, the equation (180) may be replaced 
by 

dsd/e\d/l\ , s 



dt dt\nj dt 
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where d/di is the time derivative in a local comoving Lorentz system. Now, using e = e (0) + s = s m + ,s (1) , 
V = P(o) an d n = ?i (0) + we find from Eq. (181) 

where we neglected time derivatives of first order quantities. With the help of Eqs. (170b), (170c) and (172) we find 
that the right-hand side of Eq. (182) vanishes. Hence, 



dg( ) 
dt 



0, (183) 



implying that, in zero order, the expansion takes place without generating entropy: s (0) is constant in time. Hence, 
in view of Eq. (17), which is valid for any scalar, and Eq. (183), the first order perturbation s (1) is automatically a 
gauge-invariant quantity, i.e., s (1) = s (1) , in contrast to e (1) and which had to be redefined according to Eqs. (2). 
Apparently, the entropy per particle s (1) is such a combination of e (1) and n (1) that it need not be redefined. This 
can be made explicit by noting that in the linear approximation we are considering in this article, the second law of 
thermodynamics (180) should hold for zero order and first order quantities separately. In particular, Eq. (180) implies 



T m s (1) = -L (e m - £ -^+P±!ln m ) , (184) 

n (0) \ n (0) J 

where we neglected products of differentials and first order quantities, and where we replaced de and dn by e (1) 
and n (1) , respectively. We now note that the linear combination in the right-hand side of Eq. (184) has the property 

e m - %*±*&n m = 4 _ fei+Mn* , (185) 

^(0) ^(0) 

as may immediately be verified with the help of Eqs. (2), (170b) and (170c). The right-hand side of Eq. (185) being 
gauge-invariant, the left-hand side must be gauge-invariant. This observation makes explicit the gauge invariance 
of the first order approximation to the entropy per particle, s (1) . In order to stress the gauge invariance of the 
correction s (1) to the (constant) entropy per particle, s m , we will write sfl )y rather than s (1) : 

(186) 

From Eqs. (184)-(186) we then get 

T ( o,4 = — ( 4 ~ £( ° )(1 + W) nf^ , (187) 



(0) 



where w is the quotient of zero order pressure and zero order energy density defined by Eq. (172). Notice that for 
the internal logic of our reasoning it is not essential at all to use the second law of thermodynamics. One may simply 
consider (187) as the defining expression for a certain linear combination of and n^, and replace everywhere in 
the equations below the product T^sf^ by the right-hand side of Eq. (187), without changing anything. However, 
the second law of thermodynamics yields a physical interpretation of the particular linear combination of efl } and nfl } 
that we will encounter in the equations of the next section. In fact, we rewrite Eq. (187) in the form 

g i _ £ (0) (1 + w) gi 

i (0)«(i) ~a CT (i)> 

n ( ) 

where the gauge independent, entropy related quantity afl } is given by 



The quantity ct* 1 ) occurs as the source term in the gauge-invariant evolution equations (204) below. Equation (189) 
implies (7). 
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B. Gauge-invariant pressure perturbations 

We will now derive a gauge-invariant expression for the physical pressure perturbations. To that end, we first 
calculate the time derivative of the background pressure. From Eq. (66) we have 

P(o) = Pnh (0) +p £ i m - (190) 

where p e and p n are the partial derivatives given by Eqs. (91) and (92). Multiplying both sides of this expressions by 
{ i) /0(o) and subtracting the result from Eq. (90) we get 

P w -^ (1) =p n nfi )+ p E et 1} , (191) 

0(0) 

where we have used Eqs. (2) to rewrite the right-hand side. Since p n and p £ depend on the background quantities e (0) 
and n (0) only, the right-hand side is gauge- invariant. Hence, the quantity pfl^ defined by 

&:=P W -&6 m , (192) 

0(0) 

is gauge- invariant. We thus obtain the gauge- invariant counterpart of Eq. (90) 

Pfl =Peefl ) +p n nfi ) . (193) 

We will now rewrite this equation in a slightly different form. From Eqs. (175) and (190) we obtain (3 2 — p e + 
Pn{h (0 ) /£((»)■ Using Eqs. (170b) and (170c) we find 



p=J Ps + n ;f; v (194) 

V £(o,(l + w) 

With this expression and Eqs. (189) and (194) we can express the gauge-invariant pressure (193) in terms of the 
energy density perturbation efl } and the entropy related quantity af\ rather than and the particle number density 
perturbation nfl^ 

pfl =(3 2 ef 1) + Pn afi } , (195) 

a nice result which we will not use, however. Compare this relation with the one derived by Mukhanov et al. [4], their 
equation (5.3). 



C. Gauge-invariant temperature perturbations 

Finally, we will derive an expression for the gauge- invariant temperature perturbation TJf; with the help of Eq. (A2a). 
For the time derivative of the energy density e (0) (n (0) , T (0) ) we have 



'(0) 



de \ ( de \ 

d^) T n ^ + [df) n 



T ( o). (196) 



Replacing the infinitesimal quantities in Eq. (A2a) by perturbations, we find 

Multiplying both sides of Eq. (196) by 9 (1) /9 m and subtracting the result from Eq. (197) we get 

<-&X< + (&).( r »--b4 {m 

where we have used Eqs. (2). Hence, the quantity 



T ( f ): =T (I) -^0 (1) , (199) 
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is gauge-invariant. Thus. Eq. (198) can be written as 



(i) 



de 
dn 



«7n W ' 



(200) 



implying that Tf^ can be interpreted as the gauge-invariant temperature perturbation. An equivalent form of (200) 
is given by Eq. (8). We thus have expressed the perturbation in the absolute temperature as a function of the 
perturbations in the energy density and particle number density for a given equation of state of the form e = e(n, T) 
and p = p(n, T). This equation will be used in Sec. VIII to derive an expression for the fluctuations in the background 
temperature, St, a measurable quantity. 

Finally, we give the evolution equation for the background temperature T m (t). From Eq. (196) it follows 



-3H 



£ (0) (1 + w) 



dn J j 



-(0) 



8e\ 
dTJn 



(201) 



where we have used Eqs. (170b) and (170c). This equation will be used to follow the time development of the 
background temperature once s m (t) and n {0) (t) are found from the zero order Einstein equations. 

VIII. MANIFESTLY GAUGE-INVARIANT FIRST ORDER EQUATIONS 

The five perturbation equations (173) and (174) form a set of five equations for the five unknown quantities £ (1) , 
rifi), an d Q(i)- This system of equations can be reduced in the following way. In order to arrive at non-zero 

expressions for the physical quantities ef 1} and nf[ } we have chosen in Eqs. (25) u = 9, Eq. (26), implying that 6fl } = 0, 
Eq. (27). As a consequence, we do not need the gauge dependent quantity 6 {1) . Eliminating the quantity 6 {1) from 
equations (173) with the help of Eq. (174), we arrive at the set of four first order differential equations 



e (1) + 3H(s w +p ( i,) + £ (0 )(1 + w) 



2H ( K£<1) 



^ 3i? (Dll 



0, 



h m + 3Hn {1) + n { 



2H ( KSw 



2 



0. 



+ H(2 - 3/3 2 )^ (1) + 



1 



V 2 P (1 , 



e (0 )(l + w) a 2 



= 0, 



I 3 
2 



fc(Dll = °. 



(202a) 
(202b) 
(202c) 
(202d) 



for the four quantities £ (1) , and 3 i? ( i)||- 

The system of equations (202) is now cast in a suitable form to arrive at a system of manifestly gauge-invariant 
equations for the physical quantities efl } and nfl )y since we then can immediately calculate these quantities. Indeed, 
eliminating the quantity 6 m from Eqs. (176) with the help of Eq. (174), we get 



£(1) ^ (0) - 3e (0) (1 + w) {2Hd {1) + ± A R {1) y ) 
3 R m +3k£ (0) (1 +w) 
_ 3n (0) ( K £ (1) +2H# (1) + | 3 fl(Pl|) 
n<1) % 0) +3ne m (l + w) 



(203a) 
(203b) 



The quantities efl } and nf\ are now completely determined by the system of background equations (170)-(171) and 
the first order equations (202). 



A. Evolution equations for density perturbations 



Instead of calculating efl } and fim in the way described above, we proceed by first making explicit the gauge 
invariance of the theory. To that end, we rewrite the system of the four differential equations (202) for the gauge 
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dependent variables e w , n (1) , and 3 i? ( i)|| into a system of equations for the gauge-invariant variables ef', and nfl } . 
It is, however, of convenience to use the entropy related perturbation afl } , defined by Eq. (189), rather than the 
particle number density perturbation nfn . The result is 



£(i) + ai£ ( i) + 02£(i) — a 3< J (1) , 

erf/, = 040-fi), 



(204a) 
(204b) 



where of^ is a short hand notation for the combination of ef t) and nfl } given in Eq. (189). The derivation of these 
equations is given in detail in Appendix D. The coefficients a\, . . . ,04 occurring in Eqs. (204) are given by 



"1 



a 2 



as 



"4 



K£(o) (1 



H 



K£(0) (1 



#(4-3/3 2 )4- 3 i? (0) 



1 

3H ' 
( 



2H{1 + 3/3 2 



3 R (0) + 3k£ (0) (1 



w) 



-18H 2 



2H 2 (2-3/3 2 



rtr. 



5H 2 (l + 3/3 2 )-2H^-\ 
P 

3 i? (0) +3k£ (0) (1 + u;) 



3 i? (0) + 3k£ (0) (1 + w) 



e (0) p £n (l + w) + 



2 Pn P 

3H j3 



P 2 Pr, 



■ PePr, 



n (a)Pn 




(205a) 



,(205b) 



-3H 1 



11(a)Pn 



£(o)(l 



% 0) ,(205c) 



(205d) 



where the functions iu(4) and fi(t) are given by Eqs. (172) and (175), respectively. In the derivation of the above 
results, we used Eqs. (170) and (171). The abbreviations p n and p e are given by Eq. (91). Furthermore, we used the 
abbreviations 



d 2 p 
dn 2 ' 



Pe, 



d 2 P 

de dn 



(206) 



The equations (204) contain only gauge-invariant quantities and the coefficients are scalar functions. Thus, these 
equations are manifestly gauge-invariant. 

The equations (204) are equivalent to one equation of the third order, whereas one would expect that the four first 
order equations (202) would be equivalent to one equation of the fourth order. This observation reflects the fact that 
the solutions of the first order equations are gauge dependent, while the solutions ef^ and crfij of Eq. (204) are gauge 
independent. One 'degree of freedom', say, the gauge function ip has disappeared from the scene altogether. 

The Eq. (204b) can be solved 



At 



(207) 



Inserting this expression into Eq. (204a) we obtain the final equation for the perturbation (6). 

The equations (204) constitute the main result of this article. In view of Eq. (189), they essentially are two differen- 
tial equations for the perturbations and nfl } to the energy density s m (t) and the particle number density n (0 )(t), 
respectively, for flrw universes with k = — 1,0,4-1 . They describe the evolution of the energy density perturba- 
tion e^'j and the particle number density perturbation nfl } for flrw universes filled with a fluid which is described 
by an equation of state of the form p = p(n, e), the precise form of which is left unspecified. 



B. Evolution equations for contrast functions 



In the study of the evolution of density perturbations it is of convenience to use a quantity which measures the 
perturbation to the density relative to the background densities. To that end we define the gauge-invariant contrast 
functions 5 e and 8„ by 



S e (t,x) 



£(i) 



S n (t,x) 



n (0)(t) 



(208) 
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Using these quantities, Eqs. (204) can be rewritten as (sec Appendix D 3) 

S e + b 1 S £ + b 2 S e 



1 d 



6 n ~ 

1 + w 
ZHn m p n 

c dt \ 1 + w J e m {l + w) 
where the coefficients b\ , 62 and 63 are given by 



1 + w 



(209a) 
(209b) 



_ _ 8 _ + fe + + /I 2 H(! ± 3g ) 

-H p \3ii 3 K (0) + 3k£ (0) (1 + w) / 



(210a) 



■|/te (0) (l + iu)(l + 3w) + iJ 2 (1 - 3w + 6/3 2 (2 + 3w)) 

Ke co) (l+to) \ 3n A i? 2 (l + 6u;)(l + 3/3 2 ) 



3 7? 



3k£ (0) (1 + io) 



5 i? 



(0) I 2 



! 7? 



o, -r 3k£ (0 )(1 + w) 1 



~ P 2 I - 5 3i? (o> I . ( 210b ) 



6.3 = 



-18i/ 2 



J i? (0) + 3k£ (0 ,(1 + w) 



£(0)Psn(t+W) + — -p~ P Pn + PePn + n {0} p r , 



+ Pn 




fe( 







210c) 



In Sec. X wc use the Eqs. (209) to study the evolution of small energy density perturbations and particle number 
perturbations in flrw universes. 

The entropy perturbation (187) reads, in terms of the contrast functions (208), 



£ (0)(1 + w) 



n m T m \ 1 + w 
Finally, if we define the relative temperature perturbation St by 

TJjj (t, x 



S T {t,x) := 

then the relative temperature perturbation (8) is given by 



T m (t) ' 



(211) 



(212) 



dn J 7 



n^Sn 



(213) 



We thus have found the relative temperature perturbation as a function of the relative perturbations in the energy 
density and particle number density for an equation of state of the form e = e(n, T) and p = p(n, T) (see Appendix A). 
The quantity 5t(*, x) is a measurable quantity in the cosmic background radiation. 



IX. NON-RELATIVISTIC LIMIT IN AN EXPANDING UNIVERSE 



It is well known that if the gravitational field is weak and velocities are small with respect to the velocity of light 
(v/c — > 0), the system of Einstein equations and conservation laws may reduce to the single field equation of the 
Newtonian theory of gravity, namely the Poisson equation (28). Since in the Newtonian theory the gravitational field 
is described by only one, time- independent, potential <p(x), one cannot obtain the Newtonian limit by simply taking 
the limit v/c — > 0, since in a relativistic theory the gravitational field is described, in general, by six potentials, namely 
the components hij(x) of the metric. In this article we have used the decomposition (124). Moreover, we have shown 
in Sec. IV that hy given by (125a) describes the evolution of density perturbations. By using this decomposition 
we have reduced the number of potentials to only two, namely <j)(x) and Q(x). At first sight, it would follow from 
the evolution equations (158) for scalar perturbations that a further reduction of the number of potentials by one 
is impossible. However, we have rewritten the system of equations (158) for scalar perturbations into an equivalent 
system (202). As a result, the perturbation to the metric, h\-, enters the system (202) only via the trace 



3R w\\ 



20l fe 



l\k 



+ c 



'fe I 



(214) 
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of the perturbation to the Ricci tensor (128a). This shows explicitly that density perturbations arc, in a non-flat 
flrw universe, described by two potentials 4>{x) and (,{x). Hence, in non-flat FLRW universes the perturbation 
equations (202) do not reduce to the Newtonian theory of gravity. However, for a flat flrw universe, i.e., a universe 
characterized by k = 0, implying, in view of (77), that 



3 R (0) = 0, 



the perturbation to the Ricci scalar, Eq. (214), reduces to 



4 V 2 



3R w\\ = —3 



(215) 



(216) 



where V 2 is the usual Laplace operator. Hence, in a flat flrw universe we need only one potential <j> to describe 
density perturbations. Therefore, we consider the only candidate flrw universe in which the Newtonian limit could 
be obtained, namely a flat FLRW universe. 

Usually, one considers the non-relativistic limit in the case of a static gravitational field. In this section we will show 
that our treatment of perturbations leads to the Poisson equation also for an expanding universe, i.e., for a non-static 
gravitational field. Thus, in the approach presented in this article there is no necessity to take the static limit at all. 
In fact, we will show that the set of perturbation equations (202)-(203) reduce to Eq. (28) without the assumption 
that the expansion is negligibly slow or absent. 

Let us first consider the background equations. For a flat flrw universe, the zero order equations (170) and (171) 
reduce to 



H = -±K£ (0) (1 + w), 
-3He m (l + w), 



£ (0) 

fl 



(0) — 3ii72.( ), 



and the constraint equation 



3H A 



K£r ) 



A. 



We now consider the perturbation equations. Upon substituting Eq. (216) into Eqs. (202) and putting 3 R (0) 



(217a) 
(217b) 
(217c) 

(218) 
- 0, 



we arrive at the set of perturbation equations 



e (1) + 3i7(e (1) +j> ( i)) + e«»(l + w) 



0(U 



n (1) + 3Hn {1) + n, 0) 



0(D 



t9 (1) + H(2 - 3/? 2 )tf (1) + 
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" 2H 
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2 V 2 



0, 



V 2 p (1) 



£ (0 )(1 + w) a 2 



V> + — 5-o%)(l + w)$ w = 0, 



(219a) 
(219b) 
(219c) 
(219d) 



for the four quantities £ (1) , and <p. Similarly, we obtain from Eqs. (203) 



V 2 



2H<j> 



4ttG 



(3H 2 - A)(l + ttf) 
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£ (0) (1 + w) 



( £ (i) e (i)) 



(220a) 
(220b) 



where we have eliminated with the help of Eq. (219d) and £ (0) with the help of Eq. (218). The scale factor of the 
universe a(t) follows from the Einstein equations via H := a/a. 
We now consider the limits (see Appendix C) 



U. 



(Dll 



0, p = 0. 



(221) 



In this limit, the kinetic energy is small compared to the rest energy of a particle [see, for example, Eq. (250)], so that 
the energy density of the universe in the limit (221) is 



nmc 2 , 



(222) 
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where mc 2 is the rest energy of a particle with mass m and n the particle number density. Thus, in the limit (221), 
the background equations (217)-(218) take the simple form 



H = 




(223a) 


(0) — 


—3He (0) , 


(223b) 


(0) = 


-3Hn m , 


(223c) 



while the constraint equation (218) reads 



3H 2 = ke (0) + A. (224) 



Note that in the limit (221), Eqs. (223b) and (223c) are identical, since then e m = n {0) mc 2 , in view of Eq. (222). 

We now consider the perturbation equations (219) and (220) in the limit (221). Substituting := (^jii)^ = 0, 
w := P( a) /s { a) = and p^ 1} = into these equations we find that Eq. (219c) is identically satisfied, whereas the 
remaining equations (219) reduce to 

s m + 3He m + |g (/se w -1^= 0, (225a) 

n (1 > + 3#n (1) + || (ke w - I ^) = 0, (225b) 
V 2 = 0, (225c) 



and Eqs. (220) become 




AttG 



-a'efl, (226a) 



(226b) 

Note that in the limit (221), Eqs. (225a) and (225b) are identical, since then e (1) = n (1) mc 2 , in view of Eq. (222). In 
Eqs. (225c)-(226) the quantities e (1) and n (1) do not occur. As a consequence, these equations can be solved for the 
quantities <fi, sfl } and nf',. In Appendix C we show that e (1) and n {1) arc gauge dependent also in the limit (221), so 
that Eqs. (225a) and (225b) need not be considered. 

We now consider Eqs. (225c) and (226a) in some detail. By substituting Eq. (225c) into Eq. (226a) we obtain 

W 2 4>(x) = ^a 2 (t)ef 1) (t,x). (227) 

Since V 2 is independent of time, Eq. (225c), this equation is equivalent to 



V 2 <j>(x) = ^a 2 (t p )sfi ) (t p ,x), (228) 



where t p indicates the present time. This Einstein equation can be rewritten in a form that closely resembles the 
Poisson equation, by introducing the potential tp 

*•> : = M>- <229) 

Inserting (229) into (228) we obtain the result (28): the Einstein equation (227) for the time dependent perturbation 
f?fi')(i, x) is, at fixed t = t p , identical to the usual, time independent, Poisson equation. Hence, as argued in Sec. IH, 
the perturbations sfl } and nfl } may indeed be interpreted as physical perturbations. 



X. ANALYTICAL EXAMPLES OF THE PERTURBATION THEORY 



In this section we consider two particular simple cases as an example: a flat flrw universe with a vanishing 
cosmological constant in its radiation-dominated and matter-dominated stage. In order to keep this article self- 
contained, we start with the zero order equations and their solutions, although nothing is new here. We need these 
solutions to obtain explicit forms for the first order equations. 
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A. Radiation-dominated phase of the fiat universe 

In the radiation-dominated era we have p = |e, so that, according to (172), w = |. Next to this the cosmological 
constant is put equal to zero: A = 0. Furthermore, a flat universe (k = 0) is considered, implying, with (77), that 

1. Zero-order equations and their solutions 
The zero-order equations (170) then reduce to 

H = -§K£ (0 ), (230a) 
e (0) = -4He m , (230b) 
n (0) = -3Hn w , (230c) 



while the constraint equation (171) becomes 



The general solutions of these equations are 



H 2 = |«e (0) . (231) 



H(t) = Uct)- 1 = H(t ) [ f) , (232a) 



to 



-2 



£«»(*) = ^r(ct) 2 =e ( o,(io) Q-) , (232b) 
™ ( o)0) = «(o)(*o) f4rr) ■ ( 232c ) 



a (t ) 

The initial values H(to) and e (0) (to) are related by the initial value condition (231). 
Using the definition of the Hubble parameter H := a/a we find from (232a) that 

l 

t \ 2 



a(t)=a(t )\-) , (233) 
where to is the time at which the radiation-dominated era sets in. 

2. First-order equations and their solutions 

The zero-order solutions (232) can now be substituted into the coefficients (210) of the equations (209). Recalling 
that p = in the radiation-dominated regime, p n = and p s = |, so that (3 = 1/V3, see (194). Since V 2 = V 2 for 
a flat universe, the coefficients b\, hi and &3 reduce to 

1 V 2 

h = -H, 6 a = -___ + |«e (0)) 6 3 = 0, (234) 
where we have used (231). For the first-order equations (209) this yields the simple forms 

6 E -H8 S +(-—■ + |«e (0) ^ <5 £ = 0, (235a) 
™ (*» - K) = 0, (235b) 

where H := a/a, (62). The solution of equation (235b) expresses the fact that particle number density perturbations 
are coupled to perturbations in the radiation, i.e., 

5 n {t,x) = U £ (t,x). (236) 
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This relation is, in the extreme case of radiation-domination, independent of both the nature of the particles and the 
scale of the perturbation. 

Equation (235a) may be solved by Fourier analysis of the function S e . Writing 



8 e (t,x)=5 B (t,q)d q ", 5 n (t,x)=6 n (t,q)e i «- x , 
with q = \q\ = 2tt/X, where A is the wavelength of the perturbation and i 2 = —1, yielding 

V 2 <5 e (t, x) = -q 2 5 £ (t, q), V 2 5 n (t, x) = -q 2 5 n (t, q), 
so that the evolution equation (235a) for the amplitude S e (t,q) reads 



8 e - H(t ) [-) S e + 



1 q 2 



3 a 2 (t ) V*o 



to 



<L = 0, 



(237) 
(238) 

(239) 



where we have used (231)-(232). This equation will be rewritten in such a way that the coefficients become dimcn- 
sionlcss. To that end a dimcnsionless time variable is introduced, defined by 



This definition implies 



c n dt r - 



t 

t := — , t>t . 
to 



— ) i— = [2H(t )] n i—, n=l,2,..., 
ct J dr n 1 y n dr n ' ' ' ' 



(240) 



(241) 



where we have used (232a). Using (232a), (240) and (241), equation (239) for the density contrast S s (t, q) can be 
written as 



1 „, ( H 2 



1 



4=0, 



(242) 



^4r 2r 2 

where a prime denotes differentiation with respect to r. The constant /j, t is given by 

alo) H(t )^- (243) 

The general solution of equation (242) is a linear combination of the functions J ± i (/i r -y/r)r 3 / 4 , where J + i(x) — 
y/2/ (irx) sin x and J_i(x) = y / 2/(irx) cosx are Bessel Functions of the first kind: 

5 e (r, q) = Ax(q) sin (/x r Vr) + A 2 (q) cos (firy/r) y/r, 
where the functions A\{q) and A 2 (q) are given by 



Ax(q) = 5 e (t ,q) sin fi r 
A 2 {q) = S e (t ,q) cos fi T 



COS fl v 
Mr 

sin /it r 

Mr 



S s (to,q) 
Ss(to,q) 



S £ (tp,q) 
H(t ) 

S s (t ,q) 
H(t ) 



where we have used that 



S e (t , q) = S e (t = 1, q), 6 e (t ,q) = 25' e (t = 1, q)H(t ) 



(244) 

(245a) 
(245b) 

(246) 



as follows from (241). 

For large-scale perturbations, A — ► oo, the magnitude of the wave vector \q\ = 2tt/X vanishes. Writing 
d s (t) = S E (t, q = 0) and S e (t) = 6 e (t, q = 0), we find from (243)-(245) that, for t > to, 



Se(t) 



Se(t ) 



4 fa) 

H(t ) 



t 

To 



25 £ (t 



Se(t ) 

H{to) 



i 

t \ 2 
~o 



(247) 
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The growth rates proportional to t and t 1//2 have been derived from the full set of linearized Einstein equations and 
conservation laws by a large number of researchers. See Lifshitz and Khalatnikov [2], (8.11). Adams and Canuto [17], 
(4.5b), Olson [18], page 329, Peebles [19], (86.20), and Kolb and Turner [20], (9.121). The complete solution (247), 
however, has never been presented in literature. From this solution it follows that large-scale perturbations grow only 
if the initial growth rate is large enough, i.e., 



S s (t ) > S £ (to)H(t ), 

otherwise the perturbations are decaying. 

In the small-scale limit A — > (or, equivalcntly, |q| — > oo) we find, using (243)-(245), that 



(248) 



l 



/i r - /Lt r 



1 

t \ 2 



(249) 



Due to the pressure gradients, which play a role only on small scales, a smaller growth rate than in the case of 
large-scale perturbations is found. 

In contrast to the solution which results from the 'standard equation', Eq. (273), which yields oscillating density 
perturbations (274) with a decreasing amplitude, our equation leads to a solution which shows that small-scale 
perturbations oscillate with an increasing amplitude. 

Although the growth rates t and t 1 / 2 occurring in the large-scale solution (247) are well-known, the general solution 
(244), and, in particular, the large-scale solution (247) and the small-scale solution (249), have not been found earlier. 

The approach presented in this article yields an evolution equation (235a) which differs from the standard equa- 
tion (273). This difference will be explained in detail in Sec. XI. 



B. Matter-dominated phase of the flat universe 



Once protons and electrons recombine to yield hydrogen at a temperature around 4000 K, the radiation pressure 
becomes negligible, and the equations of state reduce to those of a non-relativistic monatomic perfect gas [Weinberg 
[14], equations (15.8.20) and (15.8.21)] 

e(n, T) = nm H c 2 + §nfc B T, p(n, T) = nk B T, (250) 

where k B is Boltzmann's constant, mn the proton mass, and T the temperature of the matter. 



1. Zero-order equations and their solutions 



The maximum temperature in the matter-dominated era occurs around time to of the decoupling of matter and 
radiation: T(to) f» 4000 K. Hence, from equation (250) it follows that the pressure is negligible with respect to the 
energy density, i.e., 

1 ~ M- < 3j x 10-10, (251) 
e m-RC 1 

This implies that, to a good approximation, e (0) ±p (0) sa e (0) and £ (0) w n m m}ic 2 . Hence, in an unperturbed flat flrw 
universe, the pressure can be neglected with respect to the energy density. The above facts, yield that the Einstein 
equations and conservation laws (170)-(171) for a flat flrw universe reduce to 

H = -5«£ (0 ), (252a) 

e (0) = -3He m , (252b) 

h m = -3Hn w , (252c) 

and the constraint equation 



U 2 — 1 

n — 3^0(0), 



(253) 
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where we have put the cosmological constant A equal to zero. The general solutions of the zero-order Einstein 
equations and conservation laws are 

H(t) = l(ct)- 1 = H(t ) (J-\ , (254a) 

£(0)W = Yn {ctV2 ^ £M {£) ' (2Mb) 
n m {t) = n (0) (to) (J^j ■ (254c) 

The initial values H(to) and e {0 )(to) are related by the initial value condition (253). 

Using the definition of the Hubble parameter H := a/a, it is found from (254a) for the scale factor that 

a(t) = o(f ) (£) 3 , (255) 



where io is the time at which the matter-dominated era sets in. 



2. First-order equations and their solutions 



Using Eqs. (A5) we find from the equations of state (250) 



3 

so that, to a good approximation, using (194) 



Pe = §, _p„ = -|to h c 2 , (256) 



/?(*)^ = J|^#, (257) 
c y 3 mnr 



where v s is the speed of sound. Differentiating (257) with respect to time yields 

(258) 



P 2T (0) 

For the time development of the background temperature it is found from Eqs. (201) and (250), this time not neglecting 
the pressure with respect to the energy density, that 

t m = -2HT m . (259) 

(Note that putting w = in Eq. (201) yields the incorrect result T (0) = 0.) Combining (258) and (259) results in 

j 3 = ~H. (260) 
For the evolution of the background temperature it is found that 

T m (t)=T m (t )(^] 2 , (261) 



a(t 

where we have used that H := a/a. Using (256) it is found for equation (209b) 

~{S n -5 e ) = -2H(S n -6 e ). (262) 
c at 

A solution to this equation is 

5 n (t,x) = 5 e (t,x), (263) 
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where we have used that in the matter-dominated phase nearly all energy resides in the rest mass of the particles. 
After substituting (256), (257) and (260) into the coefficients (210) it is found that 



h = 3H, b 2 = 



| K£r (0) 



(264) 



considering that for a flat universe V 2 = V 2 . For the evolution equation for density perturbations, (209a), this results 
in the simple form 

n 2 r ,2 + 6 K£ < 



3H6 E 



0. 



(265) 



where we have used the result (263). Using (237) and (238) the evolution equation for the amplitude 5 e (t,q) can be 
rewritten as 



H 2 (t 



{ r 



having incorporated (253)-(255), (257) and (261). The constant fi m is given by 

q 1 v s (t ) 



a(t ) H(t ) c 



v s (t ) 



5 k B T (0) (t ) 
3 



?71 H 



Using the dimensionless time variable (240), it is found from (254a) that 



d n 

■ n dt n 



1 

cto 



d n 
dr" 



d" 
dr"' 



1,2, 



Using this expression, equation (266) can be rewritten in the form 



6" + -6' + 



V 2 , 



9r 8/3 9r 2 



where a prime denotes differentiation with respect to r. Hence, the general solution of equation (269) is 



S e (r,q)= B 1 (q)J + 7(2fi m T 



-1/3) 



-B 2 (<z)J_i(2// m T 



-1/3) 



-1/2 



(266) 



(267) 



(268) 



(269) 



(270) 



where B\(q) and i? 2 (<7) are arbitrary functions and J± u (x) is the Bcssel function of the first kind. In other words, in 
the matter-dominated era density perturbations oscillate with a slowly decaying amplitude. 
In the large-scale limit (\q\ — > 0), it is found that, transforming back from t to t, 



Se(t) 



5S e (to) 



24 (t ) 
H(t ) 



2 

t \ 3 



Se(t ) 



4 fa) 

H(t ) 



(271) 
(272) 



having used that 

5 s (t , q) - 5 e (t = 1, q), S £ (t , q) = |^(r = 1, q)H(t ), 

as follows from (268). The solution proportional to t 2 / 3 is a standard result. Since 8 e is gauge-invariant, the standard 

non-physical gauge mode proportional to < _1 is absent from our theory. Instead, a physical mode proportional to 
£-5/3 j g f ounc ^ 

The approach presented in this article yields an evolution equation (265) which differs from the standard equa- 
tion (279). This difference will be explained in detail in the next section. 



XI. RELATION TO EARLIER WORK 



In this section it is shown that, in contrast to the approach developed in this article, the standard theory of 
cosmological density perturbations results in decaying small-scale density perturbations in a radiation-dominated flat 
flrw universe. Furthermore, we consider the standard Newtonian treatment for density perturbations in a flat flrw 
matter-dominated universe in the light of our equations, which are refined compared to the usual equations of earlier 
treatments. 

The difference between our treatment and the standard approach is that the solutions of the equations of the 
standard approach contain the gauge function ip(x), whereas our solutions are gauge-invariant. 
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A. Radiation-dominated universe 

The standard equation for the density contrast function S which can be found e.g., in the textbook of Peacock [21], 
equation (15.25), is given by 

S + 2HS- Q^- + !«e (0) ^a = o. (273) 

This equation is derived by using special relativistic fluid mechanics and the Newtonian theory of gravity with a 
relativistic source term. In agreement with the text under equation (15.25) of this textbook, the term — |V 2 (5/a 2 
has been added. The same result, equation (273), can be found in Weinberg's classic [14], equation (15.10.57) with 
p = and v s = 1/V3- Using (237), it is found for the general solution of equation (273) that 

S(r,q) = ^±^lj 2 (^) + ip( q )ntfH(to)Y 2 (fi r ^), (274) 

Mr 

where C-y{q) and ip(q) are arbitrary functions (the integration 'constants'), r is given by (240), the constant fj, T is 
given by (243) and J u (x) and Y v {x) are Bessel functions of the first and second kind, respectively. The factors 8//x 2 
and ir [if H (to) have been inserted for convenience. The derivation of the solution (274) runs along the same lines as 
the derivation of (244). Thus, the standard equation (273) yields oscillating density perturbations with a decaying 
amplitude. 

For large-scale perturbations (|q| — > 0, or, equivalcntly, ^ r — > 0), the asymptotic expressions for the Bessel functions 
J2 and I2 are given by 

■M/W^^t, ^(jW^—t- 1 . (275) 



Substituting these expressions into (274), it is found for large-scale perturbations that 

5(t, q) = d{q)T - AH(t )^(q) T -\ (276) 

Large-scale perturbations can also be obtained from the standard equation (273) by substituting V 2 (5 = 0, i.e., 

8 + 2H5-±Ke m 5 = 0. (277) 

The general solution of this equation is, using Eqs. (232), given by (276). Thus far, the functions C\(q) and ip(q) 
are the integration 'constants' which can be determined by the initial values 6 (to, q) and 6 (to, q)- However, equation 
(277) can also be derived from the linearized Einstein equations and conservation laws for scalar perturbations (202): 
see the derivation in Appendix E. As a consequence, equation (277) is found to be also a relativistic equation and the 
quantity 5 = e^/e m is gauge-dependent. Therefore, the second term in the solution (276) is not a physical mode, 
but equal to the gauge mode 

<W(t,<?) = V>(<?) — = -AH(t )iiq)T-\ (278) 

£ (0) 

as follows from (177a), (230b) and (232a). Consequently, ip(q) should not be interpreted as an integration constant, but 
as a gauge function, which cannot be determined by imposing initial value conditions. Thus, the general solution (274) 
of the standard equation (273) depends on the gauge function %p(q) and has, as a consequence, no physical significance. 
This, in turn, implies that the standard equation (273) does not describe the evolution of density perturbations. Here 
the negative effect of the gauge function is clearly seen: up till now it was commonly accepted that small-scale 
perturbations in the radiation-dominated era of a flat flrw universe oscillate with a decaying amplitude, according 
to (274). The approach presented in this article reveals, however, that small-scale density perturbations oscillate with 
an increasing amplitude, according to (249). 

B. Matter-dominated universe 

The standard perturbation equation of the Newtonian theory of gravity is derived from approximate, non-relativistic 
equations. It reads 

/v 2 V 2 \ 

5 + 2H6-[-%— + %Ke m ) 6 = 0, (279) 
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where v s is the speed of sound. (See, e.g., Weinberg [14], Sec. 15.9, or Peacock [21], Sec. 15.2.) Using (237), the 
general solution of equation (279) is found to be 



Se(r,q) 



^- f^-D 1 (q)J +l (2f, m T 1 / 3 ) - H'(q) Vn^H(t ) J_, (2/^r 1 / 3 ) 
o y Mm 



r- 1 / 6 , (280) 



where D\(q) and ip(q) are arbitrary functions (the 'constants' of integration) and J± v (x) is the Bessel function of the 
first kind. The factors a/ 71 "/ fi^ and 4if (io)\A"At m have been inserted for convenience. The constant /i m is given 
by (267) and r is given by (240). The derivation of the solution (280) runs along the same lines as the derivation of 
(270). 

We now consider large-scale perturbations characterized by V 2 <5 = (i.e., |q| — > 0) or perturbations of all scales 
in the non-relativistic limit (i.e., v s /c — > 0). Both limits imply fi m — > 0, as follows from (267). The asymptotic 
expressions for the Bessel functions in the limit /i m — > are given by 

^(2^/3) « A j (2MnT i/») „ 3^.-./.. (281) 

15 V 7T 

Substituting these expressions into the general solution (280), results in 

S(t, q) = Di(<jr)r 2 / 3 - 3ff (t )^(g)r~ 1 . (282) 
In the limit /x m — > 0, equation (279) reduces to 

S+2H8- 5«£ (0) 5 = 0. (283) 

Using (254a)-(254c) we find the general solution (282) of this equation. Thus far, the functions Di(q) and ift(q) are 
the integration 'constants' which can be determined by the initial values 5 (to, q) and 8 (to, q). However, equation (283) 
can also be derived from the general theory of relativity, and is, as a consequence, a relativistic equation. In fact, 
this equation follows for large-scale perturbations from equations (202) (see Appendix E for a derivation), whereas 
in the non-relativistic limit it follows from equations (225a) and (225c). In both cases, however, it is based on the 
gauge-dependent quantity 8 = e (1) /e (0) . (In Appendix C it is shown that if a variable, e.g., e (1) , is gauge-dependent, 
it is also gauge-dependent in the non-relativistic limit.) As a consequence, the second term of (282) is equal to the 
gauge mode 

W(t,<?) = VX<7)— = -3H(to)iiq)T-\ (284) 

£ (0) 

as follows from (177a), (252b) and (254a). Therefore, ip(q) should not be interpreted as an integration constant, 
but as a gauge function, which cannot be determined by imposing initial value conditions. Since the solution (280) 
of equation (279) depends on the gauge function ip(q) it has no physical significance. Consequently, the standard 
equation (279) does not describe the evolution of density perturbations. Again we encounter the negative effect of 
the gauge function: up till now it was commonly accepted that for small-scale density perturbations (i.e., density 
perturbations with wave lengths much smaller than the particle horizon) a Newtonian treatment suffices and gauge 
ambiguities do not occur and that the evolution of density perturbations in the Newtonian regime is described by the 
standard equation (279). The approach presented in this article reveals, however, that in a fluid with an equation 
of state (250), density perturbations are described by the relativistic equation (265), independent of the scale of the 
perturbations. 



Appendix A: EQUATIONS OF STATE FOR THE ENERGY DENSITY AND PRESSURE 

We have used an equation of state for the pressure of the form p = p(n,e). In general, however, this equation of 
state is given in the form of two equations for the energy density e and the pressure p which contain also the absolute 
temperature T: 

e = e(n,T), p = p(n,T). (Al) 

In principle it is possible to eliminate T from the two equations (Al) to get p = p(n, e), so that our choice of the form 
p = p(n, e) is justified. In practice, however, it may in general be difficult to eliminate the temperature T from the 
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equations (Al). However, this is not necessary, since the partial derivatives p e and p n (91), the only quantities that 
are actually needed, can be found in an alternative way. From Eq. (Al) it follows 

d --=(S) T d " + (!?)„ dT - < A2a > 

*-(s) T *> + (£).*■■ < A2b » 

From Eq. (A2b) it follows that the partial derivatives (91) are 

--(*),♦(£).(£).■ 
*-(*).(£).■ 



From Eq. (A2a) it follows 



-l 



*i)A«r). • (A4 " 

(A4b) 



dT\ (de_ 
~d~e) n ~ [df 

Upon substituting the expressions (A4) into Eqs. (A3), we find for the partial derivatives defined by Eq. (91) 

where e and p are given by (Al). In order to calculate the second order derivative p nn replace p in Eq. (A5a) by p n . 
For p ee replace p in Eq. (A5b) by p e . Finally, for p en = p ne , replace p in Eq. (A5a) by p e or, equivalently, replace p 
in Eq. (A5b) by p n . 



Appendix B: GAUGE INVARIANCE OF THE FIRST ORDER EQUATIONS 

If we go over from one synchronous system of reference with coordinates x to another synchronous system of 
reference with coordinates x given by Eq. (3), we have 

£«o + W = 0, (Bl) 

as follows from the transformation rule (13) and the synchronicity condition (30). From this equation we find, using 
Eqs. (31), (33a), (33b) and (63) that ^(t,x) must be of the form 

£°=#r), C =9 lk d k ^(x) f *_ +x *( !B ) ) (B2) 

J a z {r) 

where t/)(x) and x l { x ) are arbitrary functions — of the first order — of the spatial coordinates x. The fact that the 
gauge function ip does not depend on the coordinate t anymore, as it did in general coordinates, see Eq. (5), is a 
consequence of the choice of synchronous coordinates for the original coordinates as well as for the transformed system 
of reference. It does not imply that the gauge-invariant quantities £?U and u^j, are independent of transformations 
within a plane of synchronicity only. In fact, they have been shown to be gauge-invariant, i.e., invariant under arbitrary 
infinitesimal coordinate transformations (3) , and the fact that we are considering, from now on, only transformations 
that transform a synchronous system of reference to another synchronous system of reference does not take away, of 
course, this more general property of being invariant under arbitrary transformations (3). 
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The energy density perturbation transforms according to (la), where e (0) is a solution of Eq. (170b). Similarly, the 
particle number density transforms according to (lb) where n m is a solution of Eq. (170c). Finally, as follows from 
(17), the fluid expansion scalar 9, Eq. (22c), transforms as 

§ w {x) = 9 w (x)+il>(x)0 m (t), (B3) 

where 6 (0) = 3H is a solution of Eq. (170a). 

From (17) with a = p, e, or n and (90) we find for the transformation rule for the first order perturbations to the 
pressure 

=Ps£(i) +P n n w . (B4) 
The transformation rule (18) with the four- velocity implies 

fifi) - <» " £",o, (B5) 

where we have used that vfts = S^o, Eq. (60). From the transformation rule (B5) it follows that u^ t) transforms under 
synchronicity preserving transformations (B2) as 

u° 1) (x)=u%(x)=0, (B6a) 

u\^(x) = u\ lA {x) - -J^-~g ik {x)d k ^{x). (B6b) 

We want to determine the transformation rules for and 3 i? ( i)||. Since the quantities 3 R, (46), and (51), are 
both non-scalars under general space-time transformations, the transformation rule (17) is not applicable to determine 
the transformation of their first order perturbations under infinitesimal space-time transformations — > x^, (3). 
Since u^u satisfies Eq. (158e), and since u^u transforms according to (B6), and since we know that satisfies 
Eq. (158e) with hats, one may verify, using (B4), that 

tf (1) (x):=tf (1) (*)-^M, (B7) 

satisfies Eq. (173c) with hatted quantities. The quantity is defined in analogy to in (86) 

^(i, = «>||)|fc- (B8) 

Apparently, transforms according to (B7) under arbitrary infinitesimal synchronicity preserving space-time trans- 
formations. Similarly, one may verify that 

3 i? (1) ,|(x) :=% 1}ll (x)+4H(t) ^^-I3^ (o)W(x) j ; (B9) 

satisfies Eq. (174). Apparently, (B9) is the transformation rule under arbitrary infinitesimal synchronicity preserving 
space-time transformations. An alternative way to find the results (B7) and (B9) is to write = — / and 
3 i? (1) H = — g, where / and g are unknown functions, to substitute, thereupon, and into Eqs. (173c) 

and (174), and to determine / and g such that the old equations (173c) and (174) reappear. In fact, our approach to 
define (B7), and 3 -R(i)||, (B9), is nothing but a shortcut to this procedure. 

It may now easily be verified by substitution that if e (1) , and 3 i? ( i)|| are solutions of the system 

(173)-(174), then the quantities e (1) , (la), n (1)! (lb), 8 (1) , (B3), (B7), and 3 ^ ( i)||, (B9), are, for an arbitrary 

function tp(x), also solutions of this system. 

Appendix C: GAUGE DEPENDENCE IN THE NON-RELATIVISTIC LIMIT 

If we require that the limit (221) holds true before and after a gauge transformation, we find from the transformation 
rule (B6b) for that dk^(x) = 0, or, equivalently, 

i>(x) = C, (CI) 
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Table I: The coefficients ay figuring in the equations D2. 



an =3H(l+p e ) + 
ftn (0 ) 



K£(0)(1 + W) 



211 



021 



211 



<*31 
«51 



Pi 



£(0)(i + w) a 2 

K 3 R(0) 



?,11 



■ 3 Ro 



3 i?( ) + 3k£( 0) (1 + w) 



ai2 = 3Hp n 
«22 = 3H 



032 



Vn 



"13 = £(o)(l + w) 
«23 = "(0) 
a 33 = H(2 - 3/3 2 ) 



e (0 )(l + w) a 2 
Q42 = CT43 = — 2k£ (0 )(1 + w) 

6e m H(l + w) 



a 52 = 



014 



"24 



£(0)(1 + w) 



111 



n (0) 
AH 



053 



034 = 



a 4 4 = 2H - 



3 R(a) + 3k£(o)(1 + w) 54 3 R (a) + 3k£ ( q)(1 + w) 



3 Rn 



Oil 



where C is an arbitrary, small, constant. In view of Eq. (CI), the gauge dependent functions £ (1) and n (1) transform 
under a gauge transformation in the non-relativistic limit according to [cf. Eqs. (1)] 

> £ (1 ) = £(!) + C£ (0 ), (C2a) 

n m -> n (1) = n (1) + Cra (0) . (C2b) 

Since in an expanding universe the time derivatives e m (t) ^ 0, (223b), and n m (t) 7^ 0, (223c), the functions s (1) (x) 
and n {l) {x) do still depend upon the gauge, so that these quantities have, also in the non-relativistic limit, no physical 
significance. Therefore, Eqs. (225a) and (225b) need not be considered in the non-relativistic limit: these equations 
are not part of the Newtonian theory of gravity. The only remaining physical equations are Eqs. (225c)-(226). 
Finally, it follows from equations (B2) and (CI) that the transformation (3) reduces in the non-relativistic limit to 

x° x Q - C, x l -> x l - x l (as) . (C3) 

In other words, and not surprisingly, time and space transformations are decoupled: time coordinates may be shifted, 
whereas spatial coordinates may be chosen arbitrarily. 



Appendix D: DERIVATION OF THE MANIFESTLY GAUGE-INVARIANT PERTURBATION 

EQUATIONS 

In this appendix we derive the perturbation equations (204) and the evolution equations (209). 



1. Derivation of the evolution equation for the entropy 

With the help of Eqs. (202a) and (202b) and Eqs. (170)-(172) one may verify that 

1 d ( "(o) "\ qejVi n (o)Pn \f n {0) \ 

cdt I 71 - - T^iT^H = - m I 1 " ~ T^T^H ■ (D1) 

In view of Eq. (185) one may replace £ (1) and n (1) by ef^ and nfL. Using Eq. (189) yields Eq. (204b) with coeffi- 
cient (205d) of the main text. 



2. Derivation of the evolution equation for the energy density perturbation 

We will now derive Eq. (204a). To that end, we rewrite the system (202) and Eq. (203a) in the form 

£ (1) + an£ (1) + ai 2 n m + ai3# w + «u 3 ^? ( i)|| = 0, (D2a) 

h m + a 21 s w + a 22 n {1) + a 23 i) m + a 24 3 i? ( i)|| = 0, (D2b) 

^(1) + Q!3i£(i) + a 32 n m + a 33 tf (1) + a 34 3 i? (1) || = 0, (D2c) 

3 i? (1) H + a 4 i£ ( i) + cti 2 n m + a 43 i? (1) + a 44 3 R(i)\\ = 0, (D2d) 

efi 1 ) + a 5 i£ ( i) + Oi 52 n (1) + a 53 d {1) + a 54 3 i? (1) || = 0, (D2e) 
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where the coefficients ctij (t) arc given in Tabic I. 

In calculating the coefficients 01, a 2 and 0,3, (205) in the main text, we use that the time derivative of the quotient w, 
defined by Eq. (172) is given by 



w = 3H(l + w)(w - (3 2 ), 



(D3) 



as follows from Eqs. (170b) and (175). Moreover, it is of convenience not to expand the function /3(i) defined by 
Eq. (175) since this will complicate substantially the expressions for the coefficients ai, 02 and 03. 

Step 1. We first eliminate the quantity 3 i? ( i)|| from Eqs. (D2). Differentiating Eq. (D2e) with respect to time and 
eliminating the time derivatives h m , and 3 R( 1} \\ with the help of Eqs. (D2a)-(D2d), we arrive at the equation 



efi\ +Pi£(D +P2n {1) +P3$ ( i) +P4 3 i? ( i,|| = 0, 
where the coefficients pi(t), . . . ,pi{t) are given by 

Pi = ot5i - ctbiotu - a 52 a2i - a 53 a 3i - 0(54044. 

From Eq. (D4) it follows that 

3p 1 -gi Pi P2 P3 , 



(D4) 



(D5) 



(D6) 



In this way we have expressed the quantity 3 -R ( i)|| as a linear combination of the quantities e^ 1 ,, e (1) , n (1) and 
Upon replacing 3 i? ( i)|| given by (D6), in Eqs. (D2), we arrive at the system of equations 



n w + q2if 1) + /?2i£ ( i 

3 i? (1) H +q i ef 1) +/?4i£ ( i 
£fi) + ?5efii + /?5i£(i 



+ /?42^(i 



where the coefficients qi (t) and fyj (t) are given by 



Pi 



+ /3l3^(i) 
+ #230<i) 

+ foa# w 

+ Ai30<i) 



QiPj- 



0, 

0, 
0, 
0, 
0, 



(D7a) 
(D7b) 
(D7c) 
(D7d) 
(D7e) 



(D8) 



We now have achieved that the quantity 3 -R(i)ii occurs only in Eq. (D7d). Since we are not interested in the non-physical 
quantity 3 -R(i)||, we do not need this equation any more. 

Step 2. We now proceed in the same way as in step 1: eliminating this time the quantity from the system of 
equations (D7). Differentiating Eq. (D7e) with respect to time and eliminating the time derivatives n {l) and 
with the help of Eqs. (D7a)-(D7c), we arrive at 



i^efi 1 ) + r £fi) + s i£(i) + s 2«- ( i) + s 3^(i) = 0, 
where the coefficients r(t) and Si(t) are given by 

St = $5i - 0510U - /?52/?2i - @5303i, 
r = 1 + <?5 ~ /?5l9l - /352<72 - /35393- 



(D9) 



(DlOa) 
(DlOb) 



From Eq. (D9) it follows that 
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In this way we have expressed the quantity t9 (1) as a linear combination of the quantities ef^, e 8 ',, and Upon 
replacing given by (Dll) in Eqs. (D7), we arrive at the system of equations 



£<i) 


o 95 ..gi 

- Pi 3 — £ (1) + 
S3 


(« 


-013- 

S3, 




o 15 ..gi 

- P23 — £(i) + 
S3 




r 

-^23- 
S3 


<*« 


o 95 ..gi 
— P33 — £(i) + 
S3 


(• 


r 

-P33- 
S3 



^ + (/?!! -/?13- 
S3 

/?21 ~ /323 — 
S3 



/3i3— ) = 

s 3 / 


0, 


(D12a) 


/?23— ) "(I) = 


0, 


(D12b) 


fe— J "(I) = 

S3/ 


0. 


(D12c) 


- ^43— ) 77(i) 
S3/ 


= 0, 


(D12d) 


/?53— ) 77(1, = 
S3 J 


0. 


(D12e) 



C + (An -/9 33 - 

V S3 

3 i? (1) ,l - /3 43 -4) + [?4 - (343-) eft + (/3 4 i - /3 43 - 

S3 V s 3 y v s 3 

4 - P53-4 + U (3 5 3-) ifl + (fa /3 53 - 

S3 V s 3j \ S3 

We have achieved now that the quantities and occur only in Eqs. (D12c) and (D12d), so that these equations 

will not be needed anymore. We are left, in principle, with Eqs. (D12a), (D12b) and (D12e) for the three unknown 
quantities e (1) , n (1) and efl } , but we first proceed with all five equations. 

Step 3. At first sight, the next steps would be to eliminate, successively, the quantities e (1) and 77 (1) from Eq. (D12e) 
with the help of Eqs. (D12a) and (D12b). We then would end up with a fourth order differential equation for the 
unknown quantity sf^ . 

However, it is possible to extract a second order equation for the gauge-invariant energy density from the equa- 
tions (D12). This will now be shown. Eq. (D12e) can be rewritten 

..gi .gi . gi ( , P51S3 — P53 s l \ rr>-io\ 

£?i) + aizfi) + fl2£fi) = a 3 n m + e (1) , (D13) 

V P52S3 — P53S2 / 

where the coefficients ai(t), 02(7;) and a,3(t) are given by 

01 = — 15- + (D14a) 

P53 95 

a 2 = (D14b) 

P5395 

P52S3 S 2 . . 

a 3 = . (D14c) 

P5305 95 

These are precisely the coefficients (205a)-(205c) of the main text. Furthermore, we find 

P51S3 - P53S1 _ 77 (0) 
P52S3 - /?53S2 6 {0) (l + w) 

Hence, 

. P51S3 - /?53Sl 77 (l 



(D15) 



P52S3-P53S2 £ ( o)(l + ?7j) 



(D16) 



With the help of this expression and Eq. (185) we can rewrite Eq. (D13) in the form (204a). 

The derivation of the expressions (205) from (D14) and the proof of the equality (D15) is straightforward, but 
extremely complicated. We used Maple V [22] to perform this algebraic task. 

3. Evolution equations for the contrast functions 

In this section we derive Eqs. (209). We start off with Eq. (209b). From Eq. (189) and the definitions (208) it 
follows that 

0% = 77(0, (s n - . (D17) 
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Differentiating this equation with respect to ct yields 

ato* = h {0) (5 n - ^A-) + n {0) \ 1 (t n - , (D18) 

where we have used Eq. (204b). Using Eqs. (170c), (205d) and (D17) to eliminate af\, we arrive at Eq. (209b) of the 
main text. 

Finally, we derive Eq. (209a). Upon substituting the expressions 

efl = e m 6 e , (D19a) 
efi 1 ) = e ( o)4 + £(o)4! (D19b) 
efi) = g (o)4 + 2e (0 )4 + e (0 )4, (D19c) 

into Eq. (204a), and dividing by e (0) , we find 

6i=2^.+oi, (D20a) 

£(o> 

6 2 = ^ + ai^ + a 2 , (D20b) 

£(o> £<o) 

b 3 = a 3 ^. (D20c) 

£(o> 

where we have also used Eq. (D17). These are the coefficients (210) of the main text. 

Appendix E: DERIVATION OF THE STANDARD EQUATION FOR DENSITY PERTURBATIONS 

In this appendix equations (277) and (283) of the main text are derived, for a flat flrw universe, 3 i?( j = 0, with 
vanishing cosmological constant, A = 0, using the background equations (170)-(171) and the linearized Einstein 
equations and conservation laws for scalar perturbations (202). 

From (D3) it follows that w is constant if and only if w = (3 2 for all times. Using (194) it is found for constant 
w that p n = and p e = w, i.e., the pressure does not depend on the particle number density. Consequently, in the 
derivation of equations (277) and (283) the equations (170c) for n m (t) and (202b) for n (1) (t,x) are not needed. In 
this case, the equation of state is given by 

p = we. (El) 

To derive the standard equations (277) and (283), it is required that u^n = 0, implying that 

t? (1) (t,x)=0. (E2) 

Since V 2 p ( i, = for large-scale perturbations, equation (202c) is identically satisfied. After substituting = £ (0) <5 
into equation (202a) and eliminating e (0) with the help of equation (170b), it is found that 

<5 + ^Ko,5+i 3 i? (1) ||) =0. (E3) 

Differentiating equation (E3) with respect to x° — ct and using equations (170a), (170b), (171) and (202d), yields 

5+ |(1 + w)H5 - |(1 + w) 2 K e m 6 - |(1 + w)(l - 3w) 3 R {1)1{ = 0. (E4) 

Eliminating ^ijh from equation (E4) with the help of equation (E3), yields the standard equation for large-scale 
perturbations in a flat flrw universe: 

6 + 2HS - ±K£ (0 )41 + w)(l + 3w) = 0. (E5) 

This equation has been derived by Weinberg [14], Eq. (15.10.57) and Peebles [19], Eq. (86.11). For w = | (the 
radiation-dominated era) equation (277) is found, whereas for w = (the matter-dominated era) equation (283) 
applies. 
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Using that the general solution of the background equations (170)-(171) for 3 i? (0) = 0, A = and constant w is 
given by 

£ '»«> = MlW^^C"^) < E6b » 

we find for the general solution of Eq. (E5) 

8(t, q) = E 1 (q)r^+ &w ^^+ 3 ^ - 3(1 + w)H{t )Hq)T~\ (E7) 

where E\(q) is an arbitrary function (the integration 'constant') and i()(q) is the gauge function [see (278) with w = | 
or (284) with w = 0]. The solution (E7) is exactly equal to the result found by Peebles [19], §86, Eq. (86.12). 

Finally, we note that in the derivation of (277) Peebles, §86 uses Eq. (E2) (in his notation: 9 = 0). In this case, 
Peebles' approach yields a physical mode 6 oc r and a gauge mode S oc t . For ^ Peebles finds a physical mode 
8 oc t 1 / 2 . However, in the approach presented in this article both physical modes S oc r and 6 oc t 1 / 2 [see Eq. (247)] 
follow from 07%6 second-order differential equation (235a), without taking an explicit value for t^(i)- 

Appendix F: SYMBOLS AND THEIR MEANING 



Table II: Symbols and their meaning of all quantities, except for those 
occurring in the appendices. 



Wat 1 m i~\ i~\ 1 

o y iiilhji 


a o yi i y~i it 
IVlfc! cllllll^ 


XVfcJltil tillCti H/(_| LlclLlOIl 


v 


( r)i r)o r)i \ 




V 2 / 


V • (V/) = S"f M 


(28) 


(y/r 


9 ij f\i 




V ■ v 




(142), (143) 


(V A w ( i))i 

v 2 / 




(145) 


v • (v/ = rv,,) 


(167) 


do 


derivative with respect to a: = ct 




di 


derivative with respect to x l 




hat: ~ 


computed with respect to x 


(1) 


dot: ' 


derivative with respect to a: = ct 


31) 


tilde: " 


computed with respect to three-metric §ij 


(54b), (71) 


superindex: gl 


gauge-invariant 
' contravariant derivative with respect to x k : C} k = g^i 


(2) 


superindex: '* 


C|i (125a) 


subindex: (0 ) 


background quantity 


88 


subindex: (i) 


perturbation of first order 




subindex: ( 2 > 


perturbation of second order 


(55)-(56) 


subindex: ; a 


covariant derivative with respect to x x 


(9) 


subindex: ^ 


covariant derivative with respect to x 


(36) 


subindex: , M 


derivative with respect to s M 


32 


subindex: « 


longitudinal part of a vector or tensor 


(124) 


subindex: j_ 


perpendicular part of vector or tensor 


(124) 


subindex: * 


transverse and traceless part of a tensor 


(124) 


P 


cr 1 times speed of sound: \/p(o)/e(o) 


(175) 


r A 

-L (IIS 


connection coefficients 


32) 


7 


arbitrary function 


(127a) 


5 e 


energy contrast function (gauge-invariant) 


(208) 




Kronecker delta 


(60) 


s„ 


particle number contrast function (gauge- invariant) 


(208) 


St 


temperature contrast function (gauge-invariant) 


(212) 


e 


energy density 


(22a) 




Levi-Civita tensor, ei 23 = +1 


(145) 


C,4> 


potentials due to relativistic density perturbations 


(125a), (214) 


n 


bookkeeping parameter equal to 1 


(55)-(56) 





expansion scalar in four-space 


(22c) 





expansion scalar in three-space 


(51) 


K 


8ttG/c 4 


(41) 


Xij 


time derivative metric coefficients 


(3D 
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Table II: (continued) 



~~K cosmological constant 

A wavelength 

/i thermodynamic potential 

/i m reduced wave length (matter) 

fir reduced wave length (radiation) 

£ M first order space-time translation 

7T arbitrary function 

QW £fl)/c 2 

a arbitrary scalar (12a) 

af ±) abbreviation for nf 1} - n ( o)efi)/[e ( o) (1 + w)] (188)-(189), (207) 

r dimensionless time (240) 

cf>, £ potentials due to relativistic density perturbations (125a), (214) 

ip Newtonian potential (28), (229) 

X 1 arbitrary three-vector (B2) 

ip first order time translation (gauge function) (l), (3), (5) 

uj arbitrary scalar (20) 

A a ' "® fj,... v arbitrary tensor (9) 

arbitrary rank two tensor (12c) 

a scale factor or radius of universe (54b) 

11,12,013,014 coefficients in perturbation equations (205) 

01,02,63 coefficients in perturbation equations (210) 

C arbitrary first order constant (CI) 

c speed of light 

ds 2 line element in four-space (29) 

E energy within volume V (178) 

G Newton's gravitation constant 

g n u metric tensor (13) 

p M „ metric with respect to x^ (13) 

§ij time-independent metric of three-space (54b) 

H c- x U (62) 

7i Hubble function: TL = (Aa/&t)/a 

hij minus first order perturbation of metric (82) 

J± v Bessel function of the first kind (274), (280) 

K constant relating Ricci and metric tensors (67) 

k k = -l (open), (flat), +1 (closed) (72) 

kB Boltzmann's constant (250) 

C(_ Lie derivative with respect to £ M (9) 

m rest mass of particle of cosmological fluid (222) 

rriH proton rest mass (250) 

N number of particles within volume V (178) 

iV M particle density four-flow (23) 

n particle number density (22b) 

p pressure (48) 

p e , p n partial derivatives of pressure (91) 

Pnn, Pen partial derivatives of pressure (206) 

q = \q\ magnitude of wave vector: 2tt/\ (237) 

q wave vector (237) 

3 R Ricci scalar in three-space (46) 

Ricci tensor in four-space (34) 

3 Rij Ricci tensor in three-space (37) 

5 entropy within a volume V (178) 

s entropy per particle s := S/N (179) 

T absolute temperature (178) 

T^ v energy momentum tensor (47) 

t cosmological time 

to initial cosmological time 

t p present cosmological time (28) 

U 1 * cosmological four- velocity U^U^ = c 2 (22) 

u" c" 1 ^ — 

U spatial velocity (221) 

u c~ x U ri41) 

arbitrary four- vector (12b) 

w pressure divided by energy (172) 

x space-time point x^ = (ct,x) 

X 1 * locally transformed coordinates (3) 

x spatial point x = (a; 1 , a; 2 , a; 3 ) 

Y± v Bessel function of the second kind (274) 
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Table II: (continued) 



ACKNOWLEDGMENTS 

One of the authors (pgm) is indebted to the director of the Institute of Theoretical Physics of the University of 
Amsterdam, prof. dr. ir. F. A. Bais, for moral support and travel expenses. 



[1] E. M. Lifshitz, J. Phys. X, 116 (1946). 

[2] E. M. Lifshitz and I. M. Khalatnikov, Adv. Phys. 12, 185 (1963). 

[3] L. D. Landau and E. M. Lifshitz, The Classical Theory of Fields, vol. 2 of Course of Theoretical Physics (Pergamon Press, 
Oxford, 1975), 3rd ed. 

[4] V. F. Mukhanov and H. A. Feldman and R. H. Brandenberger, Phys. Rep. 215, 203 (1992). 

[5] C.-P. Ma and E. Bertschinger, Astrophys. J. 455, 7 (1995). 

[6] L. O'Raifeartaigh and N. Straumann, Rev. Mod. Phys. 72, 1 (2000). 

[7] B. J. W. van Tent, Ph.D. thesis, Spinoza Institute, Utrecht University, The Netherlands, P. O. Box 80195, 3508 TD 

Utrecht, The Netherlands (2002). 
[8] J. Hwang and H. Noh (2005), URL http://arXiv.org/abs/astro-ph/0507159. 
[9] H. Noh and J. Hwang, Class. Quant. Grav. 22, 3181 (2005). 
[10] H. Noh and J. Hwang, Phys. Rev. D 69, 104011 (2004). 

[11] E. Bertschinger, in Cosmologie et Structure a Grande Echelle, edited by M. S. R. Schaeffer, J. Silk and J. Zinn- Justin 

(Elsevier Science, 1996), pp. 273-374, les Houches, Session LX, 1993. 
[12] J. Hwang and H. Noh (1997), URL http://arXiv.org/abs/astro-ph/9701137. 
[13] J. M. Bardeen, Phys. Rev. D22, 1882 (1980). 

[14] S. Weinberg, Gravitation and Cosmology: Principles and Applications of the General Theory of Relativity (John Wiley & 

Sons, Inc. New York, 1972), ISBN 0-471-92567-5. 
[15] J. M. Stewart, Class. Quantum Grav. 7, 1169 (1990). 

[16] J. W. York, jr., Annales de 1' Institut Henry Poincare - Section A: Physique theorique XXI, 319 (1974). 
[17] P. J. Adams and V. Canuto, Physical Review D 12, 3793 (1975). 
[18] D. W. Olson, Physical Review D 14, 327 (1976). 

[19] P. J. E. Peebles, The Large-Scale Structure of the Universe, Princeton Series in Physics (Princeton University Press, New 
Jersey, 1980). 

[20] E. W. Kolb and M. S. Turner, The Early Universe, Frontiers in Physics (Addison- Wesley, Reading, MA, 1990). 
[21] J. A. Peacock, Cosmological Physics (Cambridge University Press, 1999), ISBN 0-521-42270-1. 

[22] B. W. Char, K. O. Geddes, G. H. Gonnet, B. L. Leong, M. B. Monagan, and S. M. Watt, Maple V, Library Reference 
Manual (Springer- Verlag, Berlin, Heidelberg, New York, 1992), ISBN 3-540-97592-6. 



